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SPECIFICATION 



TITLE OF INVENTION 



Ambiguity Estimation of GNSS Signals for Three or More Carriers 
FIELD OF THE INVENTION 

The present invention relates to the field of Global Navigation Satellite Systems. More 
particularly, the present invention relates to ambiguity estimation of GNSS signals for three or 
10 more carriers. 

BACKGROUND 

Global Navigation Satellite Systems (GNSS) include the Global Positioning System (GPS), 
Glonass, and the proposed Galileo system. 

Each GPS satellite transmits continuously using two radio fi-equencies in the L-band, referred 
15 to as LI and L2, at respective fi-equencies of 1575.42 MHz and 1227.60 MHz. Two signals 
are transmitted on LI, one for civil users and the other for Department of Defense (DoD) 
authorized users. One signal is transmitted on L2, intended only for DoD authorized users. 
Each GPS signal has a carrier at the LI or L2 frequency, pseudo-random noise (PRN) codes, 
and satellite navigation data. Two different PRN codes are transmitted by each satellite: a 
20 coarse/acquisition (C/A) code and a precision (P/Y code) which is encrypted. Each C/A code 
is a unique sequence of 1023 bits, which is repeated each millisecond. 

Figure 1 schematically illustrates a typical prior-art two-carrier- frequency scenario. Receiver 
100 receives GPS signals from any number of satellites in view, such as SVl, SV2 and SVm, 
shown respectively at 1 10, 120 and 130. The signals pass through the earth's atmosphere 
25 (ionosphere and troposphere), shown schematically at 140. Each signal has two frequencies, 
LI and L2. Receiver 100 determines from the signals a respective pseudo-range, PRl, PR2, 
PRm, to each of the satellites. Atmospheric and multipath effects cause variations in the 
signal path, as indicated schematically at 150, which distort the pseudo-range determinations. 

As the C/A code can be measured with an error of about one meter, a civil receiver not using 
30 the military-only P/Y code determines position with an error in the range of meters. 

However, the phase of the LI and L2 carriers can be measured with an accuracy of 0.01 - 
0.05 cycle (2 mm-1 cm), so relative positions can be more precisely estimated with errors in 
the range of millimeters to centimeters. Techniques for precise positioning employ LI and L2 
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carrier phase and integer ambiguity resolution, a widely researched area of precise satellite 
positioning. 



Many techniques have been developed to reliably and rapidly determine the cycle ambiguity 
in carrier phase signals observed by satellite positioning equipment. Ambiguity resolution 
5 techniques normally involve the use of imambiguous code observations, and observations on 
multiple satellites to reduce the potential carrier phase ambiguities imtil a decision can be 
made to accept a single candidate. Observations of coherently-generated carrier phase signals 
further enhances the speed and reliability of ambiguity resolution. 

The process of ambiguity resolution encompasses three steps: 

10 1 . Estimation of approximate values for the ambiguity on each satellite and carrier-phase 

band, 

2. Statistical search over the potential ambiguity candidates to find an ordered list of the 
best candidates, and 

3. Validation of the top ambiguity candidate. 



15 



Obtaining good ambiguity estimates greatly reduces the effort required for statistical search 
and validation. 



The classical approach to estimation of the carrier phase ambiguities is to construct a global . 
20 filter (estimator) that includes states (unknown parameters) for: 

1 . Rover station coordinates (x,y,z), 

2. Carrier phase ambiguity terms for each satellite and each carrier frequency band, 

3. Nuisance parameters for the ionospheric bias (one per satellite), 

4. Clock and residual tropospheric bias parameters (although these states are oflen 
25 ignored when double differencing of observations is used). 

Figure 2 shows such a prior-art solution in which all ambiguities for all observed satellites 
and both fi-equencies are estimated using a single, large state vector. Raw GPS data 200 for 
LI and L2 observations of multiple satellites is prepared at 210 and supplied as a prepared 
30 data set 220 to a global Kalman fUter 230. Filter 230 supplies ambiguity estimates for the LI 
and L2 observations. 



Currently with the Global Positioning System (GPS), up to 12 satellites can be tracked above 
a user's local horizon at once. Each GPS satellite transmits on two carrier firequencies. 
35 Hence the number of states that must be updated in the filter equals, for example, 41 : 

• 3 rover coordinate states (x, y, z) 

• 12*2 ambiguity states (for dual-frequency phase observations), 

• 12 ionospheric bias parameters, 

• 1 clock and 1 tropospheric bias state. 



40 The number of floating point operations for updating an n-state Kalman filter equals 

approximately n"^ (M. Grewal et aL, Kalman filtering: theory and practice using 
MATLAB, second edition, 2001, John Wiley & Sons, New York, ISBN: 0-471-39254-5). 
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Hence, for computational speed and efficiency it is desirable to minimize the nimiber of states 
included in a single Kalman filter. 

Once the European Galileo satellite system becomes available, up to 30 satellites (nSat=30) 
may be available at once. The Galileo satellites are expected to transmit on three, or possibly 
5 four, carrier frequencies. With GPS modernization, three carrier frequencies (nFreq=3) will 
be available. See, for example, K. de Jong, Future GPS and Galileo Signals, 
GeoInformatics, September 2002 (two pages); G. Hein et aL, Galileo Frequency & Signal 
Design, GPS WORLD, June 2003, pp. 30-37; S. Cliatt, GPS Modernization, PROCEEDINGS OF 
THE GNSS 2003, April 22-25 2003, Graz Austria. If data from two reference stations 
10 (nRef=2) for three carrier frequencies are processed at once, then the number of states that 
would need to be updated in a single-filter may equal, for example: 

3 + (nSat * nFreq * nRef) + (nSat * nRef) + (3 * nRef ) = 249 states 

Even with expected enhancements to embedded computer power, it is quite possible that the 
computational burden of a single filter will be too great. For many applications, such as real- 
15 time position determination, it is critical to obtain a reliable estimate within a given time 
interval, such as once per epoch of data. Moreover, increased processing generally means 
increased power consumption in receiver equipment, which is an important consideration for 
some types of equipment such battery-powered, portable xmits. 

One approach to reducing the computational burden of a large filtering problem is to use 
20 decentralized filters. However, mathematical approaches for decentralized Kalman filtering 
as proposed by N. Carlson, Federated Square Root Filter for Decentralized Parallel 
Processing, IEEE Transactions on Aerospace and Electronic Systems, Vol.AES-26, No. 3, 
May 1990, do not address the intricacies of applying the technique to carrier-phase ambiguity 
resolution problems. 

25 Though decentralized filtering has been used in the past in control systems and estimation 
problems, and in carrier-phase ambiguity resolution for the existing two-frequency GPS 
signals, there is a need for techniques which will address future GNSS systems having three 
or more frequencies, such as Galileo and modernized GPS. 

Code and carrier phase have been used in satellite-by-satellite filters. R. Hatch proposed 
30 using L1/L2 narrow-lane code in combination with the wide-lane, dual-frequency carrier 
phase combination. See R. Hatch, The synergism of GPS code and carrier phase 

ambiguities, PROCEEDINGS OF THE 3RD INTERNATIONAL GEODETIC SYMPOSIUM ON SATELLITE 

DOPPLER Positioning, Las Cruces, New Mexico, February 1982, Vol. 2, pp 1213-1232, and 

P. MiSRA et aL, GLOBAL POSITIONING SYSTEM: SIGNALS, MEASUREMENTS, AND 

35 Performance, Ganja-Jamuna Press, 2001 , pp. 230-233. 

Minimum-error phase combination is also known, in which measurement noise is balanced 
against ionospheric bias in a minimum error sense. Detailed theory for deriving the minimum 
error phase combination is found in L. Sjoberg, The best linear combinations of LI and L2 
frequency observables in the application ofTransit/Doppler and GPSy Manuscripta 
40 Geodetica 15, 1990, pp. 17-22. 
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Figure 3 shows a prior-art solution for processing of the current two-carrier GPS signal data, 
as implemented for post-processing of two-frequency data in products of Trimble Navigation 
Limited having software package TTC 2.7. Receiver 300 supplies a GPS signal data set 305 
having observations of LI and L2 for multiple satellites. A process 310 computes coefficients 
5 315 from GPS signal data set 305. GPS signal data set 305 and coefficients 315 are supplied 
to three filter processes: a geometry filter process 320, a geometry-free ionosphere filter 
process 330, and a geometry- free and ionosphere-free code filter process 340. Geometry filter 
process 320 uses a geometry carrier-phase combination to obtain cin array 325 of ambiguity 
estimates for the geometry carrier-phase combination and associated statistical information. 

10 Ionosphere filter process 330 uses a geometry-free ionosphere carrier-phase combination to 
obtain an array 335 of ambiguity estimates for the ionosphere carrier-phase combination and 
associated statistical information. Code filter process 340 uses geometry- free and ionosphere- 
free code-carrier combinations to obtain an array 345 of ambiguity estimates for the 
geometry-free and ionosphere- free code-carrier combinations and associated statistical 

15 information. Arrays 325, 335 and 345 are supplied to a combining process 350 to obtain an 
array 355 of ambiguity estimates for all carrier-phase observations and associated statistical 
information. Array 355 is supplied to a float-solution computation process 360 to compute a 
position 365. Code filter 340 is a single, wide-lane filter which includes multipath modeling 
but does not include noise modeling. This prior-art solution is not suited to GNSS with three 

20 or more carrier frequencies, such as Galileo. 

Figure 4 schematically illustrates a three-carrier-frequency scenario, such as proposed for 
Galileo and for modernized GPS. Receiver 400 receives GNSS signals from any number of 
satellites in view, such as SVl, SV2 and SVm, shown respectively at 410, 420 and 430. The 
signals pass through the earth's atmosphere, shown schematically at 440. Each signal has 
25 three or more frequencies, fl, f2, . . .fik. Receiver 400 determines from the signals a respective 
pseudo-range, PRl, PR2, PRm, to each of the satellites. Atmospheric and multipath effects 
cause variations in the signal path, as indicated schematically at 450, which distort the 
pseudo-range determinations. 

Figure 5 shows a prior-art solution intended for the proposed Galileo system with three carrier 

30 frequencies. Receiver 500 supplies a GaUleo signal data set 505 having observations of three 
carriers for multiple satellites. A process 510 computes coefficients 515 from Gahleo signal 
data set 505. Galileo signal data set 505 and coefficients 515 are supplied to two filter 
processes: a geometry filter process 520 and a single bank of geometry- free filter processes 
530. Geometry filter process 520 uses a geometry carrier-phase combination to obtain an 

35 array 525 of ambiguity estimates for the geometry carrier-phase combination and associated 
statistical information. The single bank of geometry-free filter processes 530 employs a 
geometry-free filter per satellite, producing for all geometry-free information a single output 
array 535 of ambiguity estimates for the geometry-free combinations and associated statistical 
information. Arrays 525 and 535 are suppUed to a combining process 540 to obtain an array 

40 545 of ambiguity estimates for all carrier-phase observations and associated statistical 

information. Array 545 is supplied to a computation process 550 to compute a fixed position 
555 using integer least squares, validation, and position computation. Details of the 
geometry-free filters are found in Laboratory Experiment On Carrier Phase Positioning 
Techniques for GNSS-2 (TCAR-Test), Technical Note WP 2100: Use of Physical Space 

45 Information, ESA/ESTEC Contract No. 12.406/77/NL/DS Rider 1, a Confidential Report of 
Spectra Precision Terrasat GmbH dated December 1999. A deficiency of this proposed 
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approach to three-carrier ambiguity resolution (TCAR) is that all errors are handled as noise, 
making it computationally less efficient than desired and producing ambiguity estimates 
which are worse than desired. 

Improved methods and systems are needed for ambiguity resolution of GNSS signals having 
5 three or more carriers. 
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SUMMARY OF THE INVENTION 



Embodiments in accordance with the invention address the need for carrier-phase ambiguity 
resolution in future GNSS systems having three or more frequencies, e.g., Gahleo and 
enhanced GPS. 

5 Computationally efficient techniques are employed for obtaining good approximations of the 
carrier phase ambiguities for three or more frequency bands, greatly reducing the effort 
required for statistical search and validation is greatly reduced. 

Advantages of embodiments in accordance with the invention as compared to existing multi- 
carrier ambiguity resolution techniques include high computational efficiency for the filter 

10 components (float solution) and the ability to obtain a better knowledge of the error properties 
of individual measurements and measurement combinations. This efficiency offers one or 
more benefits. Less-performant processing components are required, enabling construction of 
receivers at lower cost and/or which consume less power. Even if the processing power is 
available, most modem CPUs allow saving power by reducing the processor clock and/or 

15 putting the processing unit in a kind of sleep mode while no computations are pending. This is 
an important cost, weight and run-time factor. And for a given processing power, 
embodiments in accordance with the invention enable use of better computational models 
than is possible with prior techniques, leading to improved system performance. 
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BRIEF DESCRIPTION OF THE DRAWING 



Figure 1 schematically illustrates a typical prior-art two-carrier-frequency scenario; 

Figure 2 shows a prior-art approach to processing of two-carrier-frequency GPS signal data; 

Figure 3 shows a further prior-art approach to processing of two-carrier-frequency GPS signal 
5 data; 

Figure 4 schematically illustrates a proposed three-carrier-frequency scenario of the prior art; 

Figure 5 shows a prior-art solution intended for the proposed Galileo system with three carrier 
frequencies; 

Figure 6 shows an architecture for GNSS position determination in accordance with an 
1 0 embodiment of the invention; 

Figure 7 shows a fiirther architecture for GNSS position determination in accordance with an 
embodiment of the invention; 

Figure 8 shows elements of of a fiiUy-factorized 3+ carrier ambiguity resolution filter in 
accordance with embodiments of the invention; 

15 Figure 9 shows a summary of the numbers of individual filters and filter banks in accordance 
with embodiments of the invention; 

Figure 10 shows structure of a geometry-free ionosphere filter bank in accordance with 
embodiments of the invention; 

Figure 1 1 A shows structure of a Quintessence filter bank in accordance with embodiments of 
20 the invention; 

Figure 1 IB shows a generalized structure in which a plurality of Quintessence filter banks is 
provided in accordance with embodiments of the invention; 

Figure 12A shows structure of a single, geometry-free code filter bank in accordance with 
embodiments of the invention; 

25 Figure 12B shows structure of a single, geometry-free and ionosphere-free code filter bank in 
accordance with embodiments of the invention; 
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Figure 13 shows structure of a number of geometry-free and ionosphere-free code fiher banks 
in accordance with embodiments of the invention; 

Figure 14 shows structure of a number of mutually-orthogonal, geometry-free and 
ionosphere-free code filter banks in accordance with embodiments of the invention; 

5 Figure 15 is a flow chart illustrating a high-level processing method in accordance with 
embodiments of the invention; 

Figure 16 is a flow chart illustrating an example of system initialization in accordance with 
embodiments of the invention; 

Figure 17 is a flow chart illustrating an example of system initialization within a single epoch 
10 in accordance with embodiments of the invention; 

Figure 18 schematically illustrates a method in accordance with embodiments of the invention 
for processing a GNSS signal data set with three or more carriers; 

Figure 19 shows an appHcation of the method of Figure 18 to obtain a computed position for a 
3+ carrier-frequency GNSS receiver; 

15 Figure 20 shows a further application of the method of Figure 18 to obtain a computed 
position for a 3+ carrier-frequency GNSS receiver; 

Figure 21 shows a further application of the method of Figure 18 to obtain computed 
positions for a 3+ carrier-frequency GNSS receiver; 

Figure 22 shows an example of a filtering architecture in accordance with embodiments of the 
20 invention; 

Figure 23 illustrates a method corresponding to the architecture of Figure 22; 

Figure 24 illustrates a further method for in accordance with embodiments of the invention; 

Figure 25 shows the architecture of a GNSS receiver in accordance with embodiments of the 
invention; 

25 Figure 26 shows the architecture of a GNSS receiver in accordance with embodiments of the 
invention; 

Figure 27 shows the architecture of a GNSS receiver in accordance with embodiments of the 
invention; 

Figure 28 shows the architecture of a dual-system GNSS receiver in accordance with 
30 embodiments of the invention; 

Figure 29 illustrates an embodiment of a reference receiver processor for preparing a 
formatted multi-band RTK data stream; 

Figure 30 illustrates an embodiment of a multi-GNSS reference receiver processor for 
preparing a formatted multi-band RTK data stream; 
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Figure 31 shows a mode of operation in accordance with embodiments of the invention in 
which a single reference station is employed; 

Figure 32 shows a mode of operation in accordance with embodiments of the invention in 
which a pluraHty of reference stations is employed in a network; and 

5 Figure 33 shows a mode of operation in accordance with embodiments of the invention in 
which a pluraUty of individual reference stations is employed. 
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DETAILED DESCRIPTION 



In the context of the present invention, the following terms are intended to have the meanings 
indicated, unless otherwise indicated by the context: 

GALILEO comprises the Galileo satellite navigation system planned for deployment by the 
5 European Space Agency. 

GLONASS comprises the Glonass satellite navigation system deployed by the Soviet Union. 

GNSS (Global Navigation SateUite System) comprises generically the GPS, GLONASS and 
GALILEO systems, similar satellite-based navigation systems which may from time to time 
become operational, and pseudolite systems. 

10 GPS comprises the NAVSTAR Global Positioning System deployed by the United States 
Department of Defense, as it currently exists and as it may exist with future modernization. 

A pseudolite, also called a pseudo-satellite, is a transmitter of GNSS-like signals. A 
pseudolite is typically terrestrial. 

A reference receiver (or reference station) is a GNSS receiver at a fixed location. 
15 A rover is a mobile GNSS receiver. 

A satellite is a transmitter of GNSS signals, and is intended to include a pseudolite. 
A user is a rover or reference receiver. 
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The following notations are used in the description below: 



Notation 


Interpretation 




Wavelength of i* carrier [m] 




Standard deviation of carrier measurement [cycles] 




Standard deviation of i* code measurement [m] 




Standard deviation of ionospheric residual [m] 




Variance of uncorrelated noise at epoch k 


CTcl 


Variance of correlated noise at epoch k 


K 


Time constant of correlated noise at epoch k 


nf . 


Nxmiber of carriers [-] 


ns 


Number of satellites tracked [-] 


ne 


Number of epochs of data [-] 


<t>i 


Carrier phase measurement for frequency i [cycles] 


Pi 


Code measurement for frequency i [m] 




Transposition of Matrix A 




Inverse of quadratic Matrix A 


V 


Column vector 




Statistics of solution sol 


I 


Identity matrix, size defined by the formula context 


\A 


Absolute value of x 


Elx\ 


Expectation value of x 


e 


Base of the natural logarithm 
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OVERVIEW 



Figure 6 is a flowchart illustrating an architecture for computing GNSS position employing 
factorized ambiguity resolution of GNSS signals for three or more carriers in accordance with 
an embodiment of the invention. A GNSS signal data set 605 is a set of observations obtained 
5 by receiving signals of multiple satellites at a receiver. GNSS signal data set 605 is supplied 
to an element 610 which prepares the data for filtering, and the resulting prepared data 615 is 
then supplied to an element 620 which applies a fuUy-factorized carrier-ambiguity-resolution 
(CAR) filter to the prepared data. 

Features and variations of the fiiUy- factorized carrier-ambiguity-resolution (CAR) filter 
10 element 620 are described in detail below. The result of applying filter element 620 to the 
prepared data is an array 625 of ambiguity estimates for all carrier-phase observations and 
associated statistical information for all transmitters (e.g. for all observed GNSS satellites 
and/or pseudolites). Array 625 is supplied to a position-computation element 630 which 
computes a receiver position 635 for the time of the observations. In the example of Figure 6, 
15 element 630 computes position 635 as a float solution. 

Figure 7 is a flowchart illustrating a further architecture for GNSS positioning employing 
factorized ambiguity resolution of GNSS signals for three or more carriers in accordance with 
an embodiment of the invention. A GNSS signal data set 705 is a set of observations obtained 
by receiving signals of multiple satellites at a receiver. Data set 705 is supphed to an element 
20 710 which prepares the data for filtering, and the resulting prepared data 715 is then supplied 
to an element 720 which applies a fuUy-factorized carrier-ambiguity-resolution (CAR) filter 
to the prepared data. 

Features and variations of the fully- factorized carrier-ambiguity-resolution (CAR) filter of 
element 720 are described in detail below. The result of applying filter element 720 to the 

25 prepared data is an array 725 of ambiguity estimates for all carrier-phase observations and 
associated statistical information for all transmitters (e.g. for all observed GNSS satellites 
and/or pseudolites). Array 725 is supplied to a position-computation element 730 which 
computes a receiver position 735 for the time of the observations. In the example of Figure 7, 
element 730 computes position 735 by applying an integer least squares process 740 to array 

30 725 and applying a validation process 750. 

Figure 8 shows the structure of a fuUy-factorized 3+ carrier ambiguity resolution filter 820 in 
accordance with embodiments of the invention, suitable for carrying out filter processes 620 
and 720. A prepared GNSS signal data set 815, such as prepared GNSS signal data set 615 or 
715, is supplied to filter 820, which contains a number of elements for carrying out sub- 

35 processes. Element 825 computes coefficients fi'om the prepared data set and passes the 
prepared data set with computed coefficients to a plurality of sub-filters. These sub-filters 
include: a single geometry filter 830; a geometry- free ionosphere filter bank 835 having one 
filter per observed satellite; one or more Quintessence filter banks 840(1) to S40(nf'2) in 
which each filter bank has one filter per observed satellite; and one or more code filter banks 

40 845(1) to 845(«/) in which each filter bank has one filter per observed satellite. Arrays 

produced by the sub-filters are supplied to a combiner 850 which provides a combined array 
855 of ambiguity estimates for all carrier phase observations with associated statistical 



Attorney Docket A-1403 



page 13 



information. Array 855, corresponding to arrays 625 and 725, is supplied to an element 860 
for position computation corresponding to processes 630 and 730. 

The number of Quintessence filter banks is two less than the nxmiber nf of carrier frequencies 
of the GNSS signal data set 815, as two combinations are used by geometry filter 830 and 
5 ionosphere-filter bank 835. For example, a single bank of Quintessence filters is provided for 
a GNSS system having three carrier fi*equencies, and two banks of Quintessence filters are 
provided for a GNSS system having four carrier frequencies. At least one combined code 
filter bank is provided; if desired, a plurality of code filter banks may be provided up to the 
number nf of carrier frequencies of the GNSS signal data set 815. For example, any number 
10 from one to three code filter banks is provided for a GNSS system having three carrier 
frequencies, and any number from one to four code filter banks is provided for a GNSS 
system having four carrier frequencies. Figure 9 shows a summary of the number of 
individual filters and filter banks in accordance with embodiments of the invention, each filter 
bank having one filter per satellite. 

15 Figure 10 shows stracture of a geometry-free ionosphere filter bank 835 in accordance with 
embodiments of the invention, for a number ns of observed satellites. Ionosphere filter bank 
835 includes a respective geometry-free ionosphere filter 1010, 1020, ... 1030, each filter 
corresponding respectively to one of the observed satellites Satl, Sat2, . . . Sat«5. Each filter is 
applied to the GNSS signal data for the corresponding satellite and supplies the result to a 

20 double-diflferencing element 1040 to obtain an array 1050 of ambiguity estimates for the 
ionosphere carrier-phase combination and associated statistical information. Array 1050 is 
supplied to a combiner, such as combiner 850, for combination with arrays from other sub- 
filters of a fixUy-factorized CAR filter such as filter 820. 

Figure 1 1 A shows structure of a Quintessence filter bank 1 100 in accordance with 
25 embodiments of the invention, for a number ns of observed satellites. Quintessence filter 

bank 1 100 includes a respective geometry- free and ionosphere-free filter 1 105, 1110, ... 

1115, each filter corresponding respectively to one of the observed satellites Satl, Sat2, . . . 

Satn5. Each filter is applied to the GNSS signal data for the corresponding satellite and 

supplies the result to a double-differencing element 1 120 to obtain an array 1 125 of ambiguity 
30 estimates for the geometry- free and ionosphere- free carrier-phase combination and associated 

statistical information. Array 1 125 is supplied to a combiner, such as combiner 850, for 

combination with arrays from other sub-filters of a fiiUy-factorized CAR filter such as filter 

820. 

In accordance with an embodiment of the invention shown in Figure 1 1 A, a single 
35 Quintessence filter bank 1 100 is provided when the GNSS system has three carrier 

frequencies. In accordance with another embodiment of the invention shown in Figure 1 IB, a 
number nf-l of Quintessence filter banks is provided where nf 'xs the number of GNSS carrier 
frequencies. For example, two Quintessence filter banks are provided in the case of four 
carrier frequencies. 

40 Figure 1 IB shows a generalized structure in which a plurality of Quintessence filter banks 
1 100, 1 130 is provide in accordance with embodiments of the invention for a number ns of 
observed satellites. Quintessence filter bank 1 100 is as shown in Figure 1 1 A. Quintessence 
filter bank 1 130 includes geometry-free and ionosphere-free filters 1 135, 1 140, ... 1 145, each 
filter corresponding respectively to one of the observed satellites Satl, Sat2, . . . Satw^. Each 
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filter is applied to the GNSS signal data for the corresponding satellite and supplies the result 
to a double-differencing element 1 150 to obtain an array 1 155 of ambiguity estimates for the 
geometry-free and ionosphere- free carrier-phase combination and associated statistical 
information. Array 1 155 is supplied to a combiner, such as combiner 850, for combination 
5 with arrays from other sub-filters of a fiiUy-factorized CAR filter such as filter 820. 

Figure 12A shows structure of a single, geometry- free code filter bank 1200 in accordance 
with embodiments of the invention, for a number ns of observed satellites. Code filter bank 
1200 includes geometry-free, but not totally ionosphere-free, filters 1205, 1210, ... 1215, 
each filter corresponding respectively to one of the observed satellites Satl, Sat2, . , . Sat«5. 
10 Each filter is applied to the GNSS signal data for the corresponding satellite and supplies the 
result to a double-differencing element 1220 to obtain an array 1225 of ambiguity estimates 
for the geometry-free and ionosphere-free code-carrier combination and associated statistical 
information. Array 1225 is supplied to a combiner, such as combiner 850, for combination 
with arrays from other sub-filters of a fiiUy-factorized CAR filter such as filter 820. 

15 Figure 12B shows structure of a single, geometry-free and ionosphere-free code filter bank 

1230 in accordance with embodiments of the invention, for a number ns of observed satellites. 
Code filter bank 1230 includes geometry-free and ionosphere-free filters 1235, 1240, ... 
1245, each filter corresponding respectively to one of the observed satellites Satl, Sat2, . . . 
Satws". Each filter is applied to the GNSS signal data for the corresponding satellite and 

20 supphes the result to a double-differencing element 1250 to obtain an array 1255 of ambiguity 
estimates for the geometry- free and ionosphere-free code-carrier combination and associated 
statistical information. Array 1225 is supplied to a combiner, such as combiner 850, for 
combination with arrays from other sub-filters of a fiiUy-factorized CAR filter such as filter 
820. 

25 Figure 13 shows structure of a number of geometry- free and ionosphere- free code filter banks 
1300, . . . 1330 in accordance with embodiments of the invention, for a number ns of observed 
satellites. The number of code filter banks 1300, . . . 1330 is any number up to the number nf 
of GNSS carrier frequencies. For example, a GNSS system having three carrier frequencies 
can have up to three code filter banks, and a GNSS system having four carrier frequencies can 

30 have up to four code filter banks. 

Referring to Figure 13, filter bank 1300 includes geometry-free and ionosphere-free filters 
1305, 1310, ... 1315, each filter corresponding respectively to one of the observed satellites 
Satl, Sat2, ... Satw^. Each filter is applied to the GNSS signal data for the corresponding 
satellite and supplies the result to a double-differencing element 1320 to obtain an array 1325 
35 of ambiguity estimates for the geometry-free and ionosphere-free code-carrier combination 
and associated statistical information. Array 1325 is supplied to a combiner, such as 
combiner 850, for combination with arrays from other sub-filters of a fiiUy-factorized CAR 
filter such as filter 820. 

Referring still to Figure 13, filter bank 1330 includes geometry-free and ionosphere-free 
40 filters 1335, 1340, . . . 1345, each filter corresponding respectively to one of the observed 
satellites Satl, Sat2, . . . Sat«5. Each filter is applied to the GNSS signal data for the 
corresponding satellite and supplies the result to a double-differencing element 1350 to obtain 
an array 1355 of ambiguity estimates for the geometry-free and ionosphere-free code-carrier 
combination and associated statistical information. Array 1355 is supplied to a combiner. 
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such as combiner 850, for combination with arrays from other sub-filters of a fully-factorized 
CAR filter such as filter 820. 

Figure 14 shows structure of a number of geometry-firee and ionosphere- free code filter banks 
1400, . . . 1430 specified as being mutually orthogonal in accordance with embodiments of the 
5 invention, for a nimiber ns of observed satellites. The number of mutually-orthogonal code 
filter banks 1400, . . . 1430 is any number up to the number nf of GNSS carrier frequencies. 
For example, a GNSS system having three carrier frequencies can have up to three code filter 
banks, and a GNSS system having four carrier frequencies can have up to four code filter 
banks. 

10 Referring to Figure 14, filter bank 1400 includes mutually-orthogonal geometry-free and 
ionosphere-free filters 1405, 1410, . . . 1415, each filter corresponding respectively to one of 
the observed satellites Satl, Sat2, . . . Sat«5. Each filter is applied to the GNSS signal data for 
the corresponding satellite and supplies the result to a double-differencing element 1420 to 
obtain an array 1425 of ambiguity estimates for the geometry- free and ionosphere-free code- 

15 carrier combination and associated statistical information. Array 1425 is supplied to a 

combiner, such as combiner 850, for combination with arrays from other sub-filters of a fiiUy- 
factorized CAR filter such as fiUer 820. 

Referring still to Figure 14, filter bank 1430 includes geometry-free and ionosphere-free 
filters 1435, 1440, . . . 1445, each filter corresponding respectively to one of the observed 

20 satellites Satl, Sat2, . . . Satw5. Each filter is applied to the GNSS signal data for the 

corresponding satellite and supplies the result to a double-differencing element 1450 to obtain 
an array 1455 of ambiguity estimates for the geometry- free and ionosphere-free code-carrier 
combination and associated statistical information. Array 1455 is supplied to a combiner, 
such as combiner 850, for combination with arrays from other sub-filters of a fixUy-factorized 

25 CAR filter such as filter 820. 

Figure 15 is a flow chart illustrjating a high-level processing method in accordance with 
embodiments of the invention. Starting at 1505, pre-processed GNSS data is received at 
1510. A fiilly-factorized 3+ carrier ambiguity resolution filter, such as fiUer 820, is initialized 
at 1515 for a first epoch and processes the data for the first epoch. Results of processing for 
30 the first epoch are provided as output at 1520. A check made at 1525 whether a pre-processed 
data set is available for a new epoch. If no, the process ends at 1530. If yes, the data set for 
the new epoch is processed at 1 535 and the results of the processing the new epoch are 
provided as output at 1520. The method repeats from 1525 as indicated. 

Figure 16 is a flow chart illustrating one of many possible examples of system initialization in 
35 accordance with embodiments of the invention. Starting at 1605, a check is made at 1610 
whether a raw GNSS data set is available for a new epoch. If no, initialization is considered 
to have failed and control returns via 1615 to restart the process at 1605. If yes, the raw 
GNSS data set is prepared at 1620 and filtered at 1625. A fiilly-factorized 3+ carrier 
ambiguity resolution combination, such as combination 820, is applied to the prepared and 
40 filtered data set at 1630. If the ambiguity estimates have not sufficiently converged, 

initialization is considered to have failed and control retums via 1615 to restart the process at 
1605. If the ambiguity estimates have sufficiently converged, initialization is considered 
successfiil and a position is computed, such as by an integer least squares calculation at 1640. 
What constitutes sufficient convergence of the ambiguity estimates to consider initialization 
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successful is a matter of design choice which may be, for example, plus or minus five or ten 
cycles of a carrier. If the ambiguity estimates are not converging, the processing resources to 
compute an integer-least-squares solution can be spared by restarting and awaiting a fresh 
data set. In any event, the ambiguity estimates will converge after a few epochs. A solution 
5 validation is performed at 1645 and a check is made at 1650 whether the validation was 
successful. If no, initialization is considered to have failed and control retums via 1615 to 
restart the process at 1605. If yes, initialization is deemed successful at 1655. 

Figure 17 is a flow chart illustrating one of many possible examples of instantaneous system 
initialization within a single epoch in accordance with embodiments of the invention. Starting 

10 at 1705, the raw GNSS data set is prepared at 1710 and filtered at 1715. A fuUy-factorized 3+ 
carrier ambiguity resolution combination, such as combination 820, is applied to the prepared 
and filtered data set at 1720. A check is made at 1725 whether initialization is possible with 
the available data. If no, initialization is considered to have failed and control retums via 
1730 to restart the process at 1705. If yes, integer ambiguities are determined, such as by an 

15 integer least squares calculation at 1735. A solution validation is performed at 1740 and a 
check is made at 1745 whether the validation was successful. If no, initialization is 
considered to have failed and control retums via 1730 to restart the process at 1705. If yes, 
initialization is deemed successful at 1750. 

Figure 1 8 schematically illustrates a method in accordance with embodiments of the invention 
20 for processing a GNSS signal data set 1805 with three or more carriers. GNSS signal data set 
1805 is optionally processed at 1810 to compute coefficients 1815 for use in the sub-filters of 
a fuUy-factorized 3+ carrier ambiguity resolution filter, such as filter 820; coefficients 1815 
may alternatively be computed in the subfilters at the cost of somewhat greater processing 
burden. The coefficients are interrelated and are orthogonal as they are computed to 
25 minimize error as discussed in more detail below. Sub-filter process 1820 applies to data set 
1805 a geometry filter such as geometry filter 830 using a geometry carrier-phase 
combination to obtain an array 1825 of ambiguity estimates for the geometry carrier-phase 
combination and associated statistical information. Sub-filter process 1830 applies to data set 
1805 a bank of ionosphere filters such as ionospher filter bank 835 using a geometry- free 
30 ionosphere carrier-phase combination to obtain an array 1835 of ambiguity estimates for the 
ionosphere carrier-phase combination and associated statistical information. 

Sub- filter process 1840 applies to data set 1805 at least one bank of Quintessence filters such 
as Quintessence filter banks 840(1). . .840(n/^2) using a geometry-free and ionosphere-free 
carrier-phase combination to obtain an array 1 845 of ambiguity estimates for the geometry- 

35 free and ionosphere-free carrier-phase combination and associated statistical information. 

The number of Quintessence filter banks is dependent on the number of carrier frequencies as 
discussed above with reference to Figures 9 and 1 IB. Sub-filter process 1850 applies to data 
set 1805 at least one code filter bank such as code filter bank 1200 or 1230 or code filter 
banks 1300 . . .1330 or 1400 . . . 1430, using a plurality of geometry- free code-carrier 

40 combinations to obtain an array 1855 of ambiguity estimates for the code-carrier 

combinations and associated statistical information for the plurality of transmitters. The 
number and characteristics of the code filter banks are described above with reference to 
Figures 9, 12A, 12B, 13 and 14. Arrays 1825, 1835, 1845 and 1855 are combined at 1860 to 
obtain a combined array 1865 of ambiguity estimates for all carrier phase observations and 

45 associated statistical information. Sub-filter processes 1820, 1830, 1840, 1850 may be carried 
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out in parallel, for example in separate threads within a processor or in separate processors, as 
desired to optimize considerations such as computational efficiency, processor power 
consvimption and/or overall processing time to position fix fi-om availability of data set 1805. 

Figure 19 shows an application of the method of Figxire 18 to obtain a computed position for a 
5 3+ carrier-fi-equency GNSS receiver 1900. GNSS signal data set 1805 is prepared by receiver 
1900 fi-om observables of 3+ carrier- frequency GNSS signals. Array 1865 is supplied to a 
process 1910 which computes a float-solution position of receiver 1900 corresponding to the 
time of acquisition of data set 1805 and provides the computed position 1920 as an output. 

Figure 20 shows a fiirther application of the method of Figure 18 to obtain a computed 
10 position for a 3+ carrier-fi-equency GNSS receiver 2000. GNSS signal data set 1805 is 

prepared by receiver 2000 from observables of 3+ earner- frequency GNSS signals. Array 
1865 is supplied to a process 2010 which determines integer ambiguities by computing an 
integer-least-squares solution. A validation process 2015 checks validity of the solution, 
which is used in a process 2020 to compute position 2025 of receiver 2000. 

15 Figure 21 shows a fiirther application of the method of Figure 18 to obtain both a float- 
solution computed position 1920 and a float-and-fixed-combination computed position 2020 
for a 3+ carrier-frequency GNSS receiver 2100. GNSS signal data set 1805 is collected by 
receiver 2100 from observables of 3+ carrier-frequency GNSS signals. Array 1865 is 
supplied to a process 1910 which computes a float-solution position of receiver 1900 

20 corresponding to the time of acquisition of data set 1805 and provides the computed position 
1920 as an output. Though unverified and potentially less precise than a validated integer- 
least-squares solution, the float-solution position is computed more quickly and can be 
supplemented with a validated integer-least-squares solution when available. Array 1865 is 
also supplied to a process 2010 which determines integer ambiguities by computing an 

25 integer-least-squares solution. A validation process 2015 checks validity of the solution. A 
process 2020 computes position from the validated solution. The computed position 2025 
supplied as an output is a validated, combination float and fixed solution. 

Decentralized estimation of carrier-phase ambiguities in accordance with embodiments of the 
invention employs the following filters: 

30 •A geometry filter, such as filter 830, which processes a minimum-error carrier-phase 

combination to estimate a single (combined) ambiguity term for each satellite, plus 
position states for the receiver; 

• A bank of ionosphere filters, such as bank 835, which estimate the ionospheric bias on 
the satellite signals using a minimum-error ionospheric carrier-phase combination; 

35 • At least one bank of Quintessence filters which process a geometry-free and 

ionosphere-free carrier-phase combination to estimate carrier-phase ambiguities 
• directly; and 

• One or more banks of code filters, up to one bank for each code band, which process 
phase combinations formed such that their ionospheric bias is equal to or near to that 

40 of the code. 



Att rney Docket A-1403 



page 18 



When multiple reference stations are processed at once, the most computationally-efficient 
approach is to implement a separate block of filters for each reference station. That is, 
geometry filter 830, ionosphere filter bank 835, one or more Quintessence filter banks 
840(1). . .840(«/-2), and one or more code filter banks 845(1), 845(2). . .845(w/) may be 
5 duplicated for each of a multiplicity of reference stations if desired for computational 
efficiency. 

The output fi-om the bank of filters is combined, in a minimum error sense, whenever 
estimates of the carrier phase ambiguities are required. Normally this would be done every 
measurement update epoch (say once per second), but if the processing platform has limited 
10 throughput, a slower update range can be used to propagate the filters and perform the 
combination. 

When multiple reference stations are involved, each reference station data stream can be 
multiplexed into the estimator to help reduce the computational load. 

A Priori Error Models 

15 A method will now be described which a priori knowledge of the error characteristics of the 
measurements used. The variance of the imcorrelated error component (noise) is employed 
together with the correlated error variance and its correlation time constant. See U. Vollath 
et aL, Network RTK Versus Single Base RTK - Understanding the Error Characteristics^ 
Proceedings of the GNSS-2002 Conference, May 2002, pp. lllA-ll^Q. 

20 This is implemented by the analysis of representative data sets. The errors are typically 
dependent on the elevation of the individual satellites and/or the carrier-to-noise ratio 
computed by the receiver tracking loops. Many methods are available in the literature. 

One method is the elevation mapping fimction approach deriving the variance cr^ at elevation 
a fi-om the variance cr^o at elevation 90°. 

25 (0-1) 

sm a 

For the time correlated errors, an exponential temporal correlation with correlation time tc is 
assumed: 

E[x{tyx{t^^t)] ^^-f 

E[x{tf]^E[x{t)f 

In the following, aul will denote the variance of an uncorrelated (white noise) error at epoch 
30 k. (Tcl is the variance of a time correlated error at epoch k with the associated time 
constant rcjt . 

Double Differences versus Single Differences 

A question of ongoing debate is the use of double differences or single differences in GNSS 
data processing. The equivalence of both approaches has been proven under the realistic 
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assumption that no precise model of the receiver clocks is available. See E. Grafarend et 
al.. Generating Classes of Equivalent Linear Models by Nuisance Parameter Elimination- 
Applications to GPS Observations, MA>aJSCRiPTA Geodetica 1 1 (1986), pp. 262-271. The 
main reason for selecting one approach or the other is the arrangement of the algorithms 
5 (software design). 

The method presented here uses single differences for the geometry- free filters, and can use 
single or double differences for the geometry filters. The ambiguity resolution step requires 
double differences for operation. 

Reference Satellite 

10 Though the ambiguity resolution process itself is not dependent on the reference satellite 
choice, a reference satellite must be selected to be able to fix ambiguities in the double 
differences. Any satellite represented in the data set can be used. Two methods of reference- 
satellite selection are common in the literature: 

1 . The highest elevation satelUte is used; 

15 2. The satellite with the lowest expected error in the measurements is selected; for 

example, the satellite with the highest signal-to-noise ratio. 

For data analysis, these reference satellite choices can provide more insights into the data 
errors. 

Data Preparation 

20 Data preparation is an optional process which pre-computes everything not dependent on the 
rover position. Operation is possible without the data preparation step at the price of 
increased processing burden in the following processes and more data to be transferred. 

To create the preprocessed carrier data for a satellite s , the frequency / and the reference 
station position jc^^ , the following formula is used: 

If a priori information about the ionosphere is also available from a model, pre-processing 
can take that into account: 

Here, Ri^^ is the geometric range between the satellite and the reference receiver position 
30 x^^ and 7?^ is the tropospheric delay between the satellite and the reference receiver position 
x^^ derived from some tropospheric model. 71^ is the tropospheric delay between the 

satellite and a good estimate (e.g., differential GNSS position) of the user position 
^user derived from the same tropospheric model. 
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Similarly, the preprocessed pseudorange measurement is: 
or, using an ionosphere model, 

p}=py'""'-pf'^-iRl -7L)+|-(^;. -^D (0-6) 

5 Li the following, all formulas refer to the preprocessed measurements. 

If no data preparation is performed, the following method descriptions have an implicit 
differencing between measurements taken at the reference receiver and the user receiver using 
the above formula. 

Geometry Filter 

10 The geometry filter, such as geometry filter 830, is a classical float solution using the 

minimum-error carrier-phase combination a^^^ (see discussion below entitled minimum-error 

geometric carrier-phase combination). The float solution can include a state for tropospheric 
residual error modeling. This is a well-known standard technique for high precision satellite 
navigation system positioning using ambiguity resolution. 

15 The float solution can be formulated in single or double differences, hi the former case, the 
double difference operator (see discussion below entitled double-difference operator) is 
applied before passing the results on to the combination operation (see discussion below 
entitled Combination). 

In embodiments of the invention, modeling of the correlated errors in the minimum-error 
20 carrier-phase combination is done by the state augmentation technique or, in the absence of 
significant noise, by the whitening-of-noise technique. 

Geometry-Free Filters 

Geometry-fi-ee filters used for the ionosphere, Quintessence and code filter banks, such as 
filter banks 835, 840(1). . .840(n/-2) and 845(1), 845(2). . .845(«/), will now be described. 
25 They implement estimating the ambiguity of a given geometry-fi-ee observable combination 
a accounting for uncorrelated noise and correlated noise with a given correlation time. 

Kalman Filter Formulation 

A geometry- free filter in accordance with embodiments of the invention is implemented as a 
Kalman filter with two states. The first state ( A^) is the ambiguity state to be estimated. The 
30 second state ( where k is the current epoch) models the time-correlated error component 
(state augmentation) with an exponential time-correlation (Gauss-Markov(l) process). 

The defining formulas for the Kalman filter (according to A. Gelb (ed.). Applied Optimal 
Estimation, The M.I.T. Press, 1992. pp. 107-113) are: 

State vector: 
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(0.7) 



State transition matrix: 



10 



<!>* = 



0 



0 ^ 



't-i 



(0.8) 



in which the upper left term 1 indicates the ambiguity is constant from epoch to epoch, k is 
the epoch number of the current epoch, k-I is the epoch number of the previous epoch, tCk is 
the time constant of correlated noise at epoch k, and the exponential term is the assvimed 
exponential temporal correlation (equation 0.2 above). 



System driving noise matrix: 



'0 
.0 



crc] 



0 ^ 

2 



(0.9) 



in which the upper left term 0 indicates that the ambiguity state (state 1) is constant 
statistically, and the lower right term crc^ is the variance of correlated noise (variance of state 
2) at epoch L 



15 Design matrix: 

H^={1 1) (0.10) 

in which the values 1 and 1 indicate the measurement contains the sum of the ambiguity state 
(state 1) and the correlated noise (state 2). 

Measurement noise matrix: 
20 Rk=(cTul) (0.11) 

is the variance of uncorrelated noise (white noise) at epoch k. 
Observation 

=<z>^ (0.12) 

is the carrier-phase measurement of epoch k for frequency/to be filtered. 
25 Filter initialization: 
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(0.13) 



in which the subscript 1 indicates starting with epoch number L 



Standard Kahnan filter algorithm (for l<k<ne): 




p-^Hl{H,^p-^Hl^R,r 
{I-K,H,)P- 



(0.14) 



5 The first two lines of 0. 14 serve to time-update the state vector and the error propagation for 
each new epoch. The last three lines of 0.14 serve to update the measurements for each new 



The first line of 0.14 produces a state vector for the current epoch A: by updating the state 
vector for the previous epoch with the state transition matrix Ok-i of the previous epoch ^-1 
10 (which in turn takes account of the time constant of the time-correlated noise). 

The second line of 0.14 produces an error propagation Pk matrix (variance-covariance matrix 
of the state) for the current epoch A: by updating the error progagation for the previous epoch 
kA with the state transition matrix O^.i of the previous epoch k-\ (which in tum takes account 
of the time constant of the time-correlated noise) and the system driving noise matrix Qk.i for 
15 the previous epoch kA (which in tum takes account of the variance of the time-correlated 



The third line of 0.14 indicates how the Kalman gain Kk for the current epoch k is related to 
the error propagation matrix Pk (which is in tum a fimction of the time-correlated noise and 
its time constant as discussed above) and to the measurement noise matrix Rk describing the 
20 variance of the uncorrelated noise. The Kalman gain Kk may be regarded as a "blending 
factor" between time updates of old information and new measurement observations. 

The fourth line of 0.14 indicates how the updated state estimate for the current epoch k is 
related to the time-updated state vector for the current epoch, the Kalman gain Kk for the 
current epoch, and the carrier-phase measurement observed for the current epoch k . 

25 The fifth line of 0. 14 is an error propagation matrix Pk^ (variance-covariance matrix of the 
state) for the ciuxent epoch k 

The Kalman filter may be implemented using the Bierman UD-Filter (see G. Bierman, 
Factorization Methods for Discrete Sequential Estimation^ Academic Press, 1977) or any other 
nimierically-stabilized implementation of the Kalman filter algorithm. 



epoch. 



noise). 
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Alternate Formulation 



In case of uncorrelated errors that are very small compared to the correlated errors, a simpler 
filter can be applied in accordance with embodiments of the invention. See the * Vhitening of 
noise" approach in G. Bierman, Factorization Methods for Discrete Sequential Estimation, 
Academic Press, 1977, and the "differencing approach" in A. Gelb, (ed.). Applied Optimal 
Estimation^ The M.I.T. Press, 1977, pp. 133-136. This applies to the ionosphere and 
Quintessence filter banks. 



Filter initialization: 



fl, =0 
9,=0 



^1 



(0.15) 



1 0 Filter algorithm {\<k<ney. 



cor^ = e 

Ct = con • - 



v^=[\-cor^^)'<Tcl 



(0.16) 



(1-c.)' 

qk=qk-x+- — — 



^k 

Xk = — 

9k 

Where Ck is a modified correlation coefficient, 4 is a modified variance (uncertainty), is the 
uncertainty in the measurement, Uk is accumulated weighted information, qk is an indicator of 
certainty with each new epoch (accumulated weights, the inverse of variance), the term a* / qk 
15 is the ambiguity, and the inverse of is the variance of ambiguity. 

Banks of Filters 

Geometry-fi-ee fihers are implemented as banks of ns filters, one filter per satellite tracked. 

The result of a geometry-fi-ee filter is determined by differencing to a specified reference 
satellite, usually highest elevation or lowest error. See discussion of reference satellite above. 
20 The choice of reference satellite does not influence the results. The results are given as 
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double differences to the chosen reference satellite. See discussion of double difference 
operator below. 



Ionosphere Filter 

The ionosphere filters, such as in ionosphere filter bank 835, are geometry- free filters using 
5 the minimum error ionospheric carrier-phase combination 5.^^^ . See discussion of minimum- 
error ionospheric carrier-phase combination below. 

Quintessence Filters 

For nf fi-equencies, /i/-2 geometry-fi-ee Quintessence filters, such as in Quintessence filter 
banks 840(1). . .840(«/^2), are implemented using the Quintessence carrier phase combinations 
10 Uq^, See discussion of Quintessence carrier-phase combinations below. 

Code Filters 

For nf fi-equencies, nf geometry-fi-ee code filters, such as in code filter banks 845(1), 
845(2). . .845(/z/), are implemented using the code-carrier combinations 5^^ . See discussion 
of code/carrier combinations below. 



15 Double Difference Operator 

Results of the geometry-free filters are given as single differences between measurements 
taken at the user (e.g., a rover) and measurements taken at a reference receiver. To convert 
the single differences into double differences, a double difference operator is appUed before 
the filter results are combined, e.g., before the results are supplied to combiner 850. A 
20 double-difference operator Mdif is for example: 



-^-1 1 0 0 
-10 10 



-1 



i : 0 

0 0 0 



• 0 0^ 

• 0 0 

•• 0 1 



(0.17) 



in which the number of columns is the number of satellites and the number of rows is one less 
than the number of satellites. In this example it is assumed that the first satellite is the 
reference satellite and that measurements for the other satellites are subtracted fi-om the 
25 measurements for the reference satellite. 

The floating-value double-difference ambiguity state estimate is the product of the double- 
difference operator M</,y-and the single-differences: 



(0.18) 



The double-difference variance-covariance matrix of the result Qdif is the product of the 
30 single-difference variance-covariance matrix Qsd and the double-difference operator and the 
transposition of the double-difference operator: 
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Qaif=M,^-Q^-Ml,, (0.19) 

Statistical values, used to validate results as discussed below, are identical for the double 
differences and the single differences: 

Xlf=x]a (0.20) 

5 Another statistical value, degrees of freedom for the geometry filter, can also be used in 
validation, e.g., to compute probability that the best solution fotmd is the correct solution: 

df,if=df^+\ (0.21) 

Similarly, degrees of freedom for geometry-free filters is given as 

dU=df^+\-ne (0.22) 

10 in which ne is the number of epochs. 
Coefficient Determination 

The following sections describe how to derive the different code and carrier combination 
coefficients for the individual filters. 

Properties of the Combinations 
15 The measurement combinations presented in the following have the properties: 

• Minimum-error geometric, minimum-error ionospheric and Quintessence 
combinations are pair-wise uncorrelated; 

• Correlation of the code combinations to the minimum-error geometric, minimum-error 
ionospheric and Quintessence combinations can be neglected as the code multipath is 

20 two to three orders of magnitude higher than the carrier-phase multipath. 

Error Models for the Combinations 

For every combination 5^^ presented, the error characteristics of uncorrelated noise crw?^ , of 
correlated noise crc?^ and time constant tc^^ can be computed using the following formulas: 



(<7(^^ ... 0 ^ 



^2 



5, (0.23) 



0 ... a-jJ^j 

25 wherein the diagonal matrix carries individual variances per frequency, and 
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a. 



tc ^-^ 

f r 



9afnf 



a. Y (a, ^ 



(0.24) 



For the code-carrier combinations, the error characteristics are practically identical to the 
error characteristics of the code measurement used. 

Minimum-Error Geometric Carrier Phase Combination 

See the general discussion of the geometry filter above. The minimum-error carrier-phase 
combination is used in the geometric float-solution filter. The carrier-phase combination used 
there is defined as: 









T 















(0.25) 



wherein (pg^o is the linear combination of carrier phases, a^eoj . . . %eo/are coefficients 
(defined below), and 0/ . . . 0«/are the original measurements at respective carrier frequencies 
1 through n . 



With the following notations: 



B =-A 



1 



(0.26) 



where Ai, X2 ■■■ ^nf are the wavelengths of the respective carrier frequencies 1 through n , 



Qgeo 



(0.27) 



where a(p2 through a(pnfZXQ the standard deviations of the measurements for carrier frequencies 2 
through w, in cycles, and 
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(0.28) 



The coefficients are defined by: 



a 



geoa 



5 for carrier fi-equencies 2 through n, and by: 



geo 



a 



geo, 2 



\J^gfiO,nf J 



(0.30) 



for carrier frequency 1 . 

Minimum-Error Ionospheric Carrier-Phase Combination 

The minimum-error ionospheric carrier-phase combination is used in the ionosphere filters. 
10 The combination is defined as: 



T'lono lono 



9\ 



a- 

\ lono^nf J 



(0.31) 



wherein (piono is the linear combination of carrier phases, aiono,i . • . atono/^e coefficients 
(defined below), and <pj ... 4>nf^^ the original measurements at respective carrier firequencies 
1 through n . 

1 5 With the following notations: 



^iono 



1 


1 ^ 















'0^ 

.1> 



(0.32) 
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10 



1 


1 ] 


-1 


r 1 



























1 



V 



(0.33) 



I 0 ait>l^ 



0 



G2 

in which the coefficients are defined by: 



(0.34) 



(0.35) 



\f^ionOtnf J 



(0.36) 



^7 ^ 



(0.37) 



Quintessence Carrier-Phase Combinations 

The Quintessence carrier-phase combinations are used in the Quintessence fihers. For 

nf carrier frequencies, k = \,-'-,nf — 2 Quintessence carrier phase combinations are defined 

by: 







^«a/ 
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(a \ 













(0.38) 



With the following notations: 



' 1 


1 
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^0^ 










0 
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1 






0 























(0.39) 
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in which the coefficients are defined by: 



"a.i 



(0.40) 



(a ^ 

"ft .3+* 




^0> 









(0.41) 



Code-Carrier Combinations 

Code-carrier combinations are used in the code filters. For every frequency k one 
combination is defined by: 



n =ar 



With the following notations: 



"n„i 



A = 



Pk 



( J_ _1^ 



r 1 A 



5„ = 



' 1 


1 > 


-1 / 





















_1_ 



1 



0n, = 



o-jZi,^ 0 ^ 
I 0 ai^lj 



0 



in which the coefficients are defined by: 

= -(^i ein, -Bn. +e2n, Gin. 



(0.42) 



(0.43) 



(0.44) 



(0.45) 



(0.46) 



(0.47) 



Attorney Docket A-1403 



page 30 



f n \ 



a, 



(0.48) 



The discussion above is directed to the general case in which measurements are made for n carrier 
frequencies. Examples are provided below for the simpler cases of 3 carrier frequencies and 4 carrier 
frequencies. 

Coefficient Example; 3 Frequencies 

The following example specifies the coefficients for the three-carrier case. 
Minimum-Error Geometric Carrier-Phase Combination 

The minimum-error carrier-phase combination is used in the geometric float solution filter. 
The carrier phase combination used there is defined as: 



—7' 



UA 






T 
















(0.49) 















With the following notations: 



^ 1 ^ 



B =-Ji 



(0.50) 



Qgeo 



^<Tll>l 0 ~ 

. 0 a4>l 



(0.51) 



A') 



(0.52) 



in which the coefficients are defined by: 

= -[<J<f^^-Ko-Ko^Qs^ +<^lKo-^l.y -(cT^' -B^^+ctI-R^^) (0.53) 



geo,2 



\^&o:3 J 



geo A T geo 

Minimum-Error Ionospheric Carrier-Phase Combination 



(0.54) 
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The minimum-error ionospheric carrier-phase combination is used in the ionosphere fiUers. 
The combination is defined as: 



T'lono iono 



iono, 3 J 



With the following notations: 



^ 1 1 
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0 0-^3^ 
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in which the coefficients are defined by: 



^ iono,2 J 



Quintessence Carrier-Phase Combination 
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(0.55) 



(0.56) 



(0.57) 



(0.58) 
(0.59) 

(0.60) 
(0.61) 



(0.62) 



the coefficients are defined by: 
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r \ 
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^0^ 
0 



(0.63) 



Code-Carrier Combinations 

Code-carrier combinations are used in the code filters. For each fi-equency ^ = 1,. . .p3 ne 
combination is defined by: 



With the following notations: 



Pk = 



a 



n.,1 



-Pk 



.A. .h. 
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r 1 N 
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1 > 




f ^ 1 
























K ^ ) 



o-j^,' 0 ^ 

L 0 C7^l^ 



10 



in which the coefficients are defined by: 

«n,.3 =-« -ein, -^n. +02^.)'' •< -gin. "^n. 



15 



Coefficient Example: 4 Frequencies 

The following example specifies the coefficients for the four-carrier case. 



(0.64) 



(0.65) 



(0.66) 

(0.67) 
(0.68) 

(0.69) 
(0.70) 
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Minimum-Error Geometric Carrier-Phase Combination 



10 



The minimum-error carrier-phase combination is used in the geometric float-solution filter. 
The carrier-phase combination used there is defined as: 



(P = a 

rgeo geo 







geo, I 
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^geo,2 










^geo,3 










\J^geoA J 







5 With the following notations: 



B 



Qgeo 



0 0 



0 o^/ 0 

L 0 0 <T^^) 



■"■geo ^ geo 



geo,3 



(0.71) 



(0.72) 



(0.73) 



(0.74) 



in which the coefficients are defined by: 

— {y<l>l 'Ko'Ko-^Qgeo ^^'Keo-KoY * ^ ' ' ^ '^,..) (0-75) 



geo 



^geo,2 
geOf^ 
\^geoA J 



(0.76) 



Minimum-Error Ionospheric Carrier-Phase Combination 

The minimum-error ionospheric carrier-phase combination is used in the ionosphere filters. 
The combination is defined as: 
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— r 



iono,\ 



a. 



iono^l 



a;. 



^3 



Vr4/ Xy-iono^ J \YAJ 



With the following notations: 



1 

















f(\\ 



lono 



a(pl 0 ^ 

0 <j<pl 



in which the coefiGcients are defined by: 



fa. ^ 

K^^iono,2 J 



(a. ,^ 

ioiio,3 



10 Quintessence Carrier-Phase Combinations 
Fork=l: 



7:T 



a, 



a. 



a ,2 



(0.77) 



(0.78) 



(0.79) 



(0.80) 



(0.81) 



(0.82) 



(0.83) 



(0.84) 



the coefficients are defined by: 
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1 1 



A 
1 



^0^ 
0 



(0.85) 



(0.86) 



For k=2: 



<3, 



e2.2 



"a ,3 



<l>y 

1^4 j 



(0.87) 



the coefficients are defined by: 
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1 
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.0, 



(0.88) 



Code-Carrier Combinations 

Code-carrier combinations are used in the code filters. For every fi-equency A: = 1, . . . , 4 one 
combination is defined by: 



a, 



*n»,3 



^^n,.4y 



Pk 



(0.89) 



With the following notations: 



4i = 



1 

















r 1 \ 



(0.90) 
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10 



15 



1 



QK = 



0 



I 0 CT^^J 
<T<I>1) 



<t4>1 0 ^ 



(0.91) 

(0.92) 
(0.93) 



the coefficients are defined by: 



n.,3 



(0.94) 



(0,95) 



Combination of Filter Output Arrays 

The combination operation, such as in combiner 850, computes the complete floating solution 
for all satellites and all frequencies from the individual filter results. The complete floating 
solution consists of the variance/covariance matrix, the float solution vector, the statistics 
and the overdetermination statistics. 

Solution Vector 

The solution vector contains the floating-value integer ambiguities for nf frequencies and ns 
satellites (giving ns-l satellite double-differences) in the following form: 

^1,1 



X = 



N 



(0.96) 



Augmented Coefficient Vectors 

The coefficient vectors are given for each individual satellite respective individual double differences. 
The coefficient vectors for the full set of satellites are defined as: 
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a = 



a. 



0 



a, 



0 

0 

0 
0 
0 



... 0 

... 0 

0 0 



0 



0 
0 
0 



0 
0 
a. 



0 a 



(0.97) 



"fj 



E.g. for three frequencies and four satellites / three double differences: 



a = 
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0> 
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0 


«3 
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«1 


0 


0 
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0 




0 


0 


0 




0 


0 


«2 


.0 


0 





(0.98) 



The matrix 5 is a «/row vector and nf •(/25'-l)x(«5-l) matrix. In the following, a denotes 
the augmented coefficient vectors . 

Variance/Covariance Matrix 

The combined variance/co variance matrix is: 



^ nf nf-2 



(0.99) 



the inverse of the sum of the estimates jfrom the geometry filter, ionospheric filter, 
10 Quintessence filter and code filter. 

Float Solution Vector 

The combined vector of floating ambiguities is: 



comb 



= Q. 



comb 



^ ^ nf-2 

Sgeo Q'Jo '^geo-^^iono'QioL '^Lo +Zi^n, ' Q'' ' ^U, + Z ^Q, ' Q'l ' 



i (0.100) 
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^ Statistics 



The statistics of the combined solution is: 



comb 





nf nf-2 



(0.101) 



The statistics provide a measure of how well the estimates fit together. 
Overdetermination produces a discrepancy between the estimates and the model. The x^ 
statistics are a weighted sum of all discrepancies between the actual measurements and the 
final model. 

Overdetermination Statistics 

The overdetermination statistics, also known as degrees-of-fireedom, is computed as: 



Integer Least Squares 

Integer Least Squares denotes a family of well-known algorithms to derive the set of integer 
ambiguities x that minimizes the function 



based on the differences between a candidate integer solution and the combined vector of 
floating ambiguities (see 0.99 above). 

This optimum set of integer ambiguities is used as the ambiguity resolution result. Validation 
is carried out as explained in the next paragraph. 

In one embodiment, the LAMBDA algorithm is used for best computational speed and efficiency. See 
P. JOOSTEN et al, GNSS Three Carrier Phase Ambiguity Resolution using the LAMBDA-method, 
Proceedings of the GNSS 1999. Other algorithms are also applicable and may be used if desired, 
such as Hatch, FARA, etc. A survey of ambiguity resolution techniques is found, for example, in P. 
MisRA et al. Global Positioning System: Signals, Measurements, and Performance, Ganja- 
Jamuna Press, 2001, Chapter 6, pp. 209-254, and in B. HOFMANN-Wellenhof et aL, GPS THEORY 
and Practice, Springer-Verlag, Fiflh Edition, 2001, pp. 213-248. 




(0.102) 



'comb 



(0.103) 



Validation 
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Validation assures that the optimum set of integer ambiguities x^^^ is the correct solution with 

a sufficient confidence. The literature gives several statistical tests that provide the probability 
of the correctness of the solution. A discussion of ambiguity validation techniques is found, 
for example, at B. Hofmann-Wellenhof et aL, GPS THEORY and Practice, Springer- Verlag, 
5 Fifth Edition, 200 1 , pp. 247-248. 

Curtin University Discrimination Test 

This test makes use of the difference between the best and the second best solution of the 
integer-least-squares problem as well as the variance/covariance matrix. See J. Wang et al., A 
discrimination test procedure for ambiguity resolution on-the-fly, JOURNAL OF GEODESY (1998) 72, 
10 pp. 644-653. 

Ratio Test / Fisher Test 

This test uses the ratio of the second best to the best solution. Either this ratio must exceed a 
predefined value (e.g. 2, or 5, or some other chosen value) or an additional statistical test 
(Fisher Test) is applied. See J. Wang et al., A discrimination test procedure for ambiguity 
15 resolution on-the-fly. Journal of Geodesy (1998) 72, pp. 644-653. 

ADOP Test 

This test is exclusively based on the variance/covariance matrix. The probability of the 
correctness of the solution is computed without inspecting the actual solution. See P.O. 
Teunissen et al.y Ambiguity dilution of precision: Definition, Properties and Application^ 
20 PROCEEDINGS OF THE ION GPS-97, 16-19 September 1997, Kansas City, USA, pp. 891-899. 

Computational Advantages 

Processing advantages for embodiments of the invention are given for very simple estimators 
(not modeling correlated errors) as well as for sophisticated estimators using state 
augmentation techniques. The Kalman filter with n states has a computational burden 
25 proportional to . A notation for this is: the n -state Kahnan filter has a computational load 
of 0{n^) 

Simple Case 

The following comparison is based on embodiments of the invention in which time-correlated 
effects are not modeled. Though these are not optimum operating modes, the comparison is 
30 given here to show a worst-case scenario for improvements offered by embodiments of the 
invention. All filters are set-up for double-differencing, which is also inferior with respect to 
the flexibility of the implementation. 

Full Filter of the Prior Art 

The prior-art fiill filter approach uses the following states: 
35 • 3 position states 

• nf • {ns - 1) ambiguity states 
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The resulting order of computational load is 

0([3 + «/.(«5-l)f) (0.104) 
Filters in Accordance with Embodiments of the Invention 

Filters and filtering methods in accordance with embodiments of the invention have the 
5 following computational burden: 

Geometry Filter 

• 3 position states 

• ns-l ambiguity states 

o([2 + «5f) (0.105) 

10 Ionosphere Filters 

If the differencing technique filter is used the ionosphere filter bank consists of ns-l separate 
1 -state Kalman filters. 

o((/w-l).[l]') (0,106) 

Quintessence Filters 

15 If the differencing technique filter is used, one Quintessence filter bank consists of 

n5-l separate 1 -state Kahnan filters and nf -2 banks of Quintessence filters are used: 

0[{nf-^2)(ns-l)[lf) (0.107) 

Code Filters 

The simple code filters use an ambiguity state only. 

20 o((«;y-l).[l]') (0.108) 
Overall Computational load 

Overall, the total computational load of the geometry filter, ionosphere filters, Quintessence 
filters and code filters for the simple case (not modeling correlated errors) is: 

o([2 + «5f +(«5-l)[lf +(w/-2).(/M-l).[lf +(«5-^^ 
25 (0.109) 

Comparison 
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The following table gives some examples of the computational load and the improvement 
factor of the simple case in accordance with embodiments of the invention over the full filter 
fomiulation of the prior art. 



No. Satellites 


Improvement 2 freq 


Improvement 3 freq 


Improvement 4 freq 


5 


3.8 


9.5 


19.0 


6 


4.2 


11.0 


22.8 


7 


4.6 


12.4 


26.1 


8 


4.8 


13.5 


28.9 


9 


5.1 


14.5 


31.4 


10 


5.3 


15.4 


33.6 



5 Complete modeling case 

The following comparison is based on embodiments of the invention in which time-correlated 
effects are modeled and using single differences. 

Full Filter 

The full filter approach uses the following states: 
10 • 3 position states 

• 1 receiver clock error state 

• nf ns ambiguity states 

• nf ns code multipath states 

• nf ns carrier multipath states 
15 • nf ns ionosphere states 

The resulting order of computational load is 

0([4 + 4.«/ «^]') (0.110) 
Filters in Accordance with Embodiments of the hivention 

Filters and filtering methods in accordance with embodiments of the invention have the 
20 following computational burden: 

Geometry Filter 
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• 3 position states 

• 1 receiver clock error state 

• ns ambiguity states 

• ns multipath states 

5 o([4 + 2-«5f) (0.111) 

Ionosphere Filters 

If the differencing technique filter is used the ionosphere filter bank consists of ns separate 1- 
state Kalman filters. 

o(«5-[lf ) (0.112) 

10 Quintessence Filters 

If the differencing technique fiher is used, one Quintessence fiher bank consists of 
ns separate 1 -state Kalman filters and nf -2 banks of Quintessence filters are used: 

o((«/-2).«5.[lf) (0.113) 

Code Filters 

15 For best results, the code filters use an ambiguity state and a multipath state to model noise 
and correlated errors simultaneously. 

o(«5[2f ) (0.114) 

Overall Computational load 

Overall, the total computational load of the geometry filter, ionosphere filter, Quintessence 
20 filters and code filters for embodiments of the invention in which time-correlated effects are 
modeled and using single differences is: 

O^^A + l ns"^ +ns\\f '^{nf'-2)'m [^^ 
Comparison 

The following table gives some examples of the load and the improvement factor of fiiU 
25 * modeling in accordance with embodiments of the invention over the prior-art fiiU filter 
formulation. 
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No. Satellites 


Improvement 2 freq 


Improvement 3 freq 


Improvement 4 freq 


5 


31 


94 


212 


6 


34 


106 


240 


7 


37 


115 


264 


8 


39 


124 


284 


9 


41 


131 


302 


10 


43 


137 


317 



Mixed Comparison 

Processing in accordance with embodiments of the invention offer an advantage for three or 
5 more carrier frequencies even when comparing fuU-modeUng embodiments in accordance 
with the invention with the simple-modeUng case for the prior-art fiill fiUers. That is, 
embodiments in accordance with the invention can offer improved processing (full-modeling 
vs. simple-modeling) while presenting a reduced computational burden over the prior-art full 
filters for three or more carrier frequencies. 

10 



No. Satellites 


Improvement 2 freq 


Improvement 3 freq 


Improvement 4 freq 


5 


0.5 


1.2 


2.5 


6 


0.5 


1.4 


2.9 


7 


0.6 


1.6 


3.3 


8 


0.6 


1.7 


3.7 


9 


0.6 


1.8 


4.0 


10 


0.7 


1.9 


4.3 



For two frequencies, the transition to full modeling has only a small impact when carried out 
using techniques in accordance with embodiments of the invention as compared with the 
simple-modeling full filter approach of the prior art. For three and four frequencies, there is 
even a speed improvement. 



Attorney Docket A-1403 



page 44 



Multiple GNSS Processing 

The methods presented are applicable to data from multiple GNSS. For example, a GPS data 
set and a Galileo data set may be processed in accordance with embodiments of the invention. 
The frequencies and therefore the combinations are determined separately for each satellite 
5 navigation system. In addition, a reference satellite is chosen for each GNSS, having a double 
difference operator application for each system. The performance of a combined satellite 
navigation system is significantly improved compared to the individual GNSS. This provides 
for availability, reliability and accuracy of the positions computed. 

Figure 22 show^s an example of a filtering architecture in accordance with embodiments of the 
10 invention. A GNSS signal data set 2205 for three or more carrier frequencies from a first 
GNSS (e.g., a Galileo data set or a modernized GPS data set) is obtained from a suitable 
GNSS receiver. A GNSS signal data set 2210 for two or more carrier frequencies from a 
second GNSS (e.g., a two-frequency or three-frequency GPS data set, a two-frequency 
Glonass data set, or a three- or four frequency Galileo data set) is obtained from a suitable 
15 GNSS receiver. An optional coefficient computation element 2215 computes coefficients for 
the filters of the first GNSS. 

GNSS data set 2205 and coefficients from element 2215, if provided, are supplied to a 
geometry-free ionosphere filter 2235, to one or more banks 2240 of Quintessence filters 
having one filter per tracked satellite of the first GNSS, and one to «/banks of code filters 

20 2245. The nxmiber of Quintessence filter banks 2240 is two less than the number of carrier 
frequencies of the first GNSS. GNSS data set 2210 and coefficients from element 2220, if 
provided, are supplied to a geometry-free ionosphere filter 2250, to a bank 2255 of 
Quintessence filters having one filter per tracked satellite of the first GNSS, and one to nf 
banks of code filters 2245. The number of Quintessence filter banks 2255 is two less than the 

25 number of carrier frequencies of the second GNSS; if the second GNSS has only two carrier 
frequencies, no Quintessence filter bank 2255 is provided. 

In the embodiment of Figure 22, a single geometry filter 2265 is common to the first GNSS 
and second GNSS as they determine the same antenna position. That is, GNSS data set 2205 
and coefficients from element 2215, if provided, are supplied to geometry filter 2265, and 
30 GNSS data set 2210 and coefficients from element 2220, if provided, are supplied to 
geometry filter 2265. 

Figure 23 illustrates a method corresponding to the architecture of Figure 22. At 2305 
geometry filter 2265 is applied to data set 2205 of the first GNSS and to data set 2210 of the 
second GNSS using a geometry carrier-phase combination to obtain an array 2308 of 

35 ambiguity estimates for the geometry carrier-phase combination and associated statistical 
information. At 2310 bank 2240 of ionosphere filters is applied to data set 2205 of the first 
GNSS using a geometry- free ionosphere carrier-phase combination to obtain an array 2315 of 
ambiguity estimates for the ionosphere carrier-phase combination and associated statistical 
information. At 2320 bank or banks 2240 of Quintessence filters are applied to data set 2205 

40 of the the first GNSS using a geometry-free and ionosphere-free carrier-phase combination to 
obtain an array 2325 of ambiguity estimates for the geometry-free and ionosphere- free 
carrier-phase combination and associated statistical information. At 2330 code filter or filters 
2245 are applied to data set 2205 of the first GNSS using a plurality of geometry-free and 
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ionosphere-free code-carrier combinations to obtain an array 2335 of ambiguity estimates for 
the code-carrier combinations and associated statistical information for the plurality of 
transmitters. 

At 2340 bank 2250 of ionosphere filters is applied to data set 2210 of the second GNSS using 
5 a geometry-free ionosphere carrier-phase combination to obtain an array 2345 of ambiguity 
estimates for the ionosphere carrier-phase combination and associated statistical information. 
At 2350 bank or banks 2255 of Quintessence filters, if provided, are applied to data set 2210 
of the the second GNSS using a geometry-free and ionosphere-free carrier-phase combination 
to obtain an array 2355 of ambiguity estimates for the geometry-free and ionosphere-free 
10 carrier-phase combination and associated statistical information. At 2360 code filter or filters 
2260 are applied to data set 2210 of the second GNSS using a plurality of geometry-free and 
ionosphere-free code-carrier combinations to obtain an array 2365 of ambiguity estimates for 
the code-carrier combinations and associated statistical information for the plurality of 
transmitters. 

15 The collection of arrays is combined at 2370 by element 2270 to produce a combined array 
2375 of ambiguity estimates for all carrier phase observations of the first GNSS and the 
second GNSS and associated statistical information. For computational efficiency, 
coefficients 2380 may be computed at 2385 for data set 2205 and coefficients 2390 may be 
computed at 2390 for data set 2210. If separately computed as shown, coefficients 2380 and 

20 2390 are supplied to the filters as indicated. 

While the architecture of Figure 22 employs a single geometry filter 2265 common to the first 
GNSS and the second GNSS, separate geometry filters may be employed for the first GNSS 
and the second GNSS if desired. Figure 24 illustrates a method in which separate geometry 
filters are employed. Items in Figures 22, 23 and 24 with like reference numerals indicate like 

25 data and processes. As shown in Figure 24, geometry filter 2410 is applied to data set 2205 of 
the first GNSS using a geometry carrier-phase combination to obtain an array 2415 of 
ambiguity estimates for the geometry carrier-phase combination and associated statistical 
information. Similarly, geometry filter 2420 is applied to data set 2205 of the second GNSS 
using a geometry carrier-phase combination to obtain an array 2425 of ambiguity estimates 

30 for the geometry carrier-phase combination and associated statistical information. Arrays 
2315, 2325, 2335, 2345, 2355, 2365, 2415 and 2425 are combined at 2430 to produce a 
combined array 2440 of ambiguity estimates for all carrier phase observations of the first 
GNSS and the second GNSS and associated statistical information. 

Mixed Numbers of Frequencies 

35 Embodiments in accordance with the invention can be applied to GNSS consisting of 

satellites transmitting different numbers of frequencies, e.g. several satellites transmit two, the 
rest three frequencies. This can be due to the capabilities of the satellites or due to limitations 
from the current transmission/reception situation. 

To process mixed numbers of frequencies, the coefficients are prepared for the lowest number 
40 of frequencies available. The additional information from the satellites transmitting more 

frequencies is taken into account by the Quintessence-filter combinations, which are available 
for at least three frequencies and are not dependent on the total number of frequencies. Also, 
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the reference satellite is chosen from the satellites with the maximum number of frequencies 
available. 

A motivation for this is that the Galileo System will provide three or four frequencies from 
the very beginning, while GPS modernization, leading eventually to three carriers (not four) 
will be introduced over many years during which some satellites will be upgraded to three 
frequencies while others will have only two frequencies available. Using two GPS carriers 
and three or four Galileo carriers will be common in the near ftiture, while using three GPS 
and three or four Galileo frequencies is a future option. 

Static Processing 

Methods and apparatus in accordance with embodiments of the invention are applicable for 
static positioning. For static positioning, the geometric float solution models one position 
coordinate set instead of one position per epoch. 

Networked GNSS Positioning 

For processing data using multiple reference stations, the ambiguities are extended to model 
all double differences from the user receiver to every reference receiver. 

Multi-Reference Station Processing 

The computational efficiency of the approach presented opens up another important 
innovation - multi-reference-station processing. RTK products on the market today only 
process a data stream from a single reference station at once. Enabling the processing of 
more than one reference station data stream simultaneously in accordance with embodiments 
of the invention offer one or more advantages: 

1 . Accuracy - error sources that are local to a single reference station tend to average out 
when multiple reference stations are processed at once, 

2. Integrity - the results obtained from a single reference station can be cross-checked 
via results from a different reference station. Today, the switching from one reference 
station to another is manual. 

3. Improved ambiguity resolution performance - the reliability and speed with which 
carrier phase ambiguities can be resolved is improved when spatial error sources are 
well modeled. If multiple reference stations are used to process GPS observations at a 
rover station, the atmospheric error sources can be better modeled than with a single 
reference. 

Multiple Rover Processing 

Processing for multiple rovers is implemented adding all ambiguities to all rovers. This 
allows embodiments of the invention to be employed in multi-rover applications such as 
determination of blade orientation of an earth-moving machine, e.g., in which a rover is 
mounted at each end of the blade. 

Precise Point Positioning 

Precise point positioning uses raw data respective differences between satellites for one 
receiver only. 
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Network Correction Processing 



Network correction processing uses ambiguity resolution for a multitude of reference stations. 
The ambiguities between all stations are resolved. This can be done efficiently using 
embodiments of the invention. All positions are pre-determined and do not have to be 
5 modeled; as the positions of the reference stations are all known, the float solution does not 
contain position states. 

Post-Processing 

Although embodiments presented are defined for real-time operation, alternate embodiments 
in accordance with the invention can also be applied by collecting the data in real-time, then 
10 processing all data collected in batch (post-processing). This applies to all processing modes 
described above. 

Instantaneous Ambignitv Resolution 

Embodiments of the invention described above are defined for multiple epochs of data used in 
the filters. Still, embodiment of the invention work for only one epoch of data 
15 (instantaneous). If the validation is successful, the precise positions can directly be computed 
from one epoch of data. This approach has the advantage of not being sensitive to cycle slips 
or loss of lock when tracking the satellites. The modified flow chart for instantaneous 
ambiguity resolution is described above with reference to Figure 17. 

Receiver Architecture 

20 Figure 25 shows the architecture of a GNSS receiver 2500 in accordance with an embodiment 
of the invention. An antenna 2505 receives GNSS signals having more than two carrier 
frequencies, such as Galileo signals having three or four carrier frequencies or modemized 
GPS signals having three carrier frequencies. The received signal passes through a low-noise 
amplifier 2510 and the amplified signal is passed by a signal splitter 2515 to a plurality of 

25 radio-frequency-band downconverters 2520, 2525, . . . 2530. The downconverters are 
supplied with a frequency reference from a reference oscillator 2535. 

Downconverter 2520 reduces a first band of the GNSS signal to a lower-frequency signal 
which is digitized by an analog-to-digital converter 2540. Downconverter 2525 reduces a 
second band of the GNSS signal to a lower-frequency signal which is digitized by an analog- 
30 to-digital converter 2545. Dovraconverter 2530 reduces an nth band of the GNSS signal to a 
lower-frequency signal which is digitized by an analog-to-digital converter 2550. 

The digitized GNSS signal of the first band is processed by a bank of detection and tracking 
channels 2560 having one channel per GNSS satellite or pseudolite tracked by receiver 2500. 
The resulting carrier-phase and pseudorange measurements for the first band are supplied to 

35 one or more processors 2575. The digitized GNSS signal of the second band is processed by 
a bank of detection and tracking chaimels 2565 having one channel per GNSS satellite or 
pseudolite tracked by receiver 2500. The resulting carrier-phase and pseudorange 
measurements for the second band are supplied to one or more processors 2575. The 
digitized GNSS signal of the nth band is processed by a bank of detection and tracking 

40 chaimels 2570 having one channel per GNSS satellite or pseudolite tracked by receiver 2500. 
The resulting carrier-phase and pseudorange measurements for the nth band are supplied to 
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one or more processors 2575. A respective downconverter, analog-to-digital converter and 
bank of detection and tracking channels is provided for each carrier-frequency band of the 
GNSS signals; for convenience of illustration only three of each are shown in Figure 25. 

Within the one or more processors 2575 is a factorized 3+ carrier filter architecture 2580 
5 which, in operation, runs a process 2585. Process 2585 processes the GNSS signal data 

comprising carrier-phase and pseudorange measurements for each of the three or more GNSS 
signal bands for each of the transmitters being tracked by receiver 2500. Filter architecture 
2580 comprises a factorized filter architecture in accordance with one or more embodiments 
of the present invention, as described for example with reference to Figures 8-14. Process 
10 2585 comprises a factorized processing method in accordance with one or more embodiments 
of the invention as described for example with reference to Figures 6-7 and 15-21. The one or 
more processors 2575 may be implemented as a single processor or as multiple processors, 
and may be programmed general-purpose processors and/or special-purpose processors. 

Process 2585 results in one or more types of output data such as: a combined array 2590 of 
15 ambiguity estimates for all carrier phase observations and associated statistical information; a 
position 2595; and/or a formatted multi-band real-time-kinematic (RTK) data stream 2598. 
Combined array 2590 corresponds, for example, to one of arrays 625, 725, 855, 1865 as 
described above. Position 2595 corresponds, for example, to one of positions 635, 735, 1920, 
2020 as described above. Data stream 2598 may be in any format suitable for use in RTK 
20 positioning equipment. 

Figure 26 shows the architecture of a GNSS receiver 2600 in accordance with an embodiment 
of the invention. An antenna 2605 receives GNSS signals having more than two carrier 
frequencies, such as Galileo signals having three or four carrier frequencies or modernized 
GPS signals having three carrier frequencies. The received signal passes through a low-noise 
25 amphfier 2610 and the amplified signal is passed to a module 2615 which performs 
downconversion from radio frequency and digitization of the multi-band signal. 

The digitized multi-band signal is processed by a bank of detection and tracking channels 
2620 having one channel per GNSS satellite or pseudolite tracked by receiver 2600. The 
resulting carrier-phase and pseudorange measurements for each band are supplied to one or 
30 more processors 2625. Reference oscillator 2622 provides a frequency reference for 
downconversion and tracking. 

Within the one or more processors 2625 is a factorized 3+ carrier filter architecture 2630 
which, in operation, runs a process 2635. Process 2635 processes the GNSS signal data 
comprising carrier-phase and pseudorange measxwements for each of the three or more GNSS 

35 signal bands for each of the transmitters being tracked by receiver 2600. Filter architecture 
2630 comprises a factorized filter architecture in accordance with one or more embodiments 
of the present invention, as described for example with reference to Figures 8-14. Process 
2635 comprises a factorized processing method in accordance with one or more embodiments 
of the invention as described for example with reference to Figures 6-7 and 1 5-2 1 . The one or 

40 more processors 2625 may be implemented as a single processor or as multiple processors, 
and may be programmed general-purpose processors and/or special-purpose processors. 

Process 2635 results in one or more types of output data such as: a combined array 2640 of 
ambiguity estimates for all carrier phase observations and associated statistical information; a 
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position 2645; and/or a formatted multi-band real-time-kinematic (RTK) data stream 2650. 
Combined array 2640 corresponds, for example, to one of arrays 625, 725, 855, 1865 as 
described above. Position 2645 corresponds, for example, to one of positions 635, 735, 1920, 
2020 as described above. Data stream 2650 may be in any format suitable for use in RTK 
5 positioning equipment. 

Figure 27 shows the architecture of a GNSS receiver 2700 in accordance with an embodiment 
of the invention. An antenna 2705 receives GNSS signals having more than two carrier 
frequencies, such as Galileo signals having three or four carrier frequencies or modernized 
GPS signals having three carrier frequencies. The received signal passes through a low-noise 
10 amplifier 2710 and the amplified signal is passed to a module 2715 which directly digitizes 
the multi-band signal. 

The digitized signals are processed by a bank of detection and tracking channels 2730 having 
one channel per GNSS satellite or pseudolite tracked by receiver 2700. The resulting carrier- 
phase and pseudorange measurements for each band are supplied to one or more processors 
15 2740. Reference oscillator 2720 provides a frequency reference for downconversion and 
tracking. 

Within the one or more processors 2740 is a factorized 3+ carrier filter architecture 2745 
which, in operation, runs a process 2750. Process 2750 processes the GNSS signal data 
comprising carrier-phase and pseudorange measurements for each of the three or more GNSS 

20 signal bands for each of the transmitters being tracked by receiver 2700. Filter architecture 
2745 comprises a factorized filter architecture in accordance with one or more embodiments 
of the present invention, as described for example with reference to Figures 8-14. Process 
2750 comprises a factorized processing method in accordance with one or more embodiments 
of the invention as described for example with reference to Figures 6-7 and 15-21. The one or 

25 more processors 2740 may be implemented as a single processor or as multiple processors, 
and may be programmed general-purpose processors and/or special-purpose processors. 

Process 2750 results in one or more types of output data such as: a combined array 2750 of 
ambiguity estimates for all carrier phase observations and associated statistical information; a 
position 2755; and/or a formatted multi-band real-time-kinematic (RTK) data stream 2760. 
30 Combined array 2750 corresponds, for example, to one of arrays 625, 725, 855, 1865 as 

described above. Position 2755 corresponds, for example, to one of positions 635, 735, 1920, 
2020 as described above. Data stream 2760 may be in any format suitable for use in RTK 
positioning equipment. 

Figure 28 shows the architecture of a GNSS receiver 2800 in accordance with an embodiment 
35 of the invention. Receiver 2800 is similar in structure to that of Figure 25, with additional 
channels for receiving and processing signals of a second GNSS. An antenna 2805 receives 
signals of a first GNSS having more than two carrier frequencies, such as Galileo signals 
having three or four carrier frequencies or modernized GPS signals having three carrier 
frequencies, and signals of a second GNSS having at least two carrier frequencies, such as 
40 GPS, modemized GPS or Glonass signals. The received signals pass through a low-noise 
amplifier 2510 and the amplified signals are passed by a signal splitter 2815 to a plurality of 
radio-frequency-band dovmconverters 2520, 2525, ... 2530 for the first GNSS and to a 
plurality of radio-frequency-band downcoverters 2820, 2825, . . .2830 for the second GNSS. 
The downconverters are supplied with a frequency reference from reference oscillators 2835. 
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Downconverter 2520 reduces a first band of the first GNSS signal to a lower-firequency signal 
which is digitized by an analog-to-digital converter 2540. Downconverter 2525 reduces a 
second band of the first GNSS signal to a lower-fi"equency signal which is digitized by an 
analog-to-digital converter 2545. Downconverter 2530 reduces an nth band of the first GNSS 

5 signal to a lower-fi-equency signal which is digitized by an analog-to-digital converter 2550. 
Downconverter 2820 reduces a first band of the second GNSS signal to a lower-fi-equency 
signal which is digitized by an analog-to-digital converter 2840. Downconverter 2825 
reduces a second band of the second GNSS signal to a lower-firequency signal which is 
digitized by an analog-to-digital converter 2845. Downconverter 2830 reduces an nth band of 

10 the second GNSS signal to a lower-fi-equency signal which is digitized by an analog-to-digital 
converter 2850, provided the second GNSS signal has more than two bands. 

The digitized signal of the first band of the first GNSS signal is processed by a bank of 
detection and tracking channels 2560 having one channel per GNSS satellite or pseudolite 
tracked by receiver 2800. The resulting carrier-phase and pseudorange measurements for the 

15 first band of the first GNSS signal are supplied to one or more processors 2875. The digitized 
signal of the second band of the first GNSS signal is processed by a bank of detection and 
tracking channels 2565 having one channel per GNSS satellite or pseudolite tracked by 
receiver 2800. The resulting carrier-phase and pseudoraiige measurements for the second 
band of the first GNSS signal are supplied to one or more processors 2875. The digitized 

20 signal of the nth band of the first GNSS is processed by a bank of detection and tracking 

channels 2570 having one channel per GNSS satellite or pseudolite tracked by receiver 2800. 
The resulting carrier-phase and pseudorange measurements for the wth band are supplied to 
one or more processors 2875. A respective downconverter, analog-to-digital converter and 
bank of detection and tracking channels is provided for each carrier-firequency band of the 

25 signals of the first GNSS; for convenience of illustration only three of each are shown in 
Figure 28. 

The digitized signal of the first band of the second GNSS signal is processed by a bank of 
detection and tracking channels 2860 having one channel per GNSS satellite or pseudolite 
tracked by receiver 2800. The resulting carrier-phase and pseudorange measurements for the 

30 first band of the second GNSS signal are supplied to one or more processors 2875. The 
digitized signal of the second band of the second GNSS signal is processed by a bank of 
detection and tracking channels 2865 having one channel per GNSS satellite or pseudolite 
tracked by receiver 2800. The resulting carrier-phase and pseudorange measurements for the 
second band of the second GNSS signal are supplied to one or more processors 2875. The 

35 digitized signal of the nth band of the second GNSS signal is processed by a bank of detection 
and tracking channels 2870 having one channel per GNSS satellite or pseudolite tracked by 
receiver 2800. The resulting carrier-phase and pseudorange measurements for the nth band 
are supplied to one or more processors 2875. A respective downconverter, analog-to-digital 
converter and bank of detection and tracking channels is provided for each carrier-frequency 

40 band of the signals of the second GNSS; though three of each are shown in Figure 28 only 
two are provided if the second GNSS has only two carrier bands. 

Within the one or more processors 2575 is a factorized 3+ carrier filter architecture 2880 
which, in operation, runs a process 2885. Process 2585 processes the signal data comprising 
carrier-phase and pseudorange measurements for each of the three or more first-GNSS-signal 
45 bands and for each of the two or more second-GNSS-signal bands for each of the transmitters 
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being tracked by receiver 2500. Filter architecture 2880 comprises a factorized filter 
architecture in accordance with one or more embodiments of the present invention, as 
described for example with reference to Figure 22. Process 2885 comprises a factorized 
processing method in accordance with one or more embodiments of the invention as 
5 described for example with reference to Figures 23 and 24. The one or more processors 2875 
may be implemented as a single processor or as multiple processors, and may be programmed 
general-purpose processors and/or special-purpose processors. 

Process 2885 results in one or more types of output data such as: a combined array 2890 of 
ambiguity estimates for all carrier phase observations and associated statistical information; a 
10 position 2895; and/or a formatted multi-band real-time-kinematic (RTK) data stream 2898. 
Combined array 2890 corresponds, for example, to one of arrays 2375 or 2440 as described 
above. Data stream 2898 may be in any format suitable for use in RTK positioning 
equipment. 

RTK Applications 

15 Embodiments of methods and apparatus in accordance with the invention may be used to 
advantage in RTK applications. 

Figure 29 illustrates an embodiment of a reference receiver processor 2905 for preparing a 
formatted multi-band RTK data stream, e.g., as indicated in Figures 25-28. Carrier-phase and 
pseudorange measurements for each of GNSS bands 1, 2, ... n are collected, e.g., in memory 
20 blocks 2910, 2915, ... 2920. The reference receiver clock supplies a time signal to a block 
2925 which provides time time tags for each epoch of measurement data, e.g., in intervals of 
0.1 second. A data set per epoch is supplied by memory blocks 2910, 2915, ... 2920 to an 
optional data compression element 2935. 

Also supplied to data compression element 2935 is a known or calculated position 2930 of the 
25 reference receiver, which is used in RTK positioning to calculate a vector firom the reference 
receiver to the rover. Data compression enables transmission overhead (time and/or cost) to 
be minimized, or enables the data rate to be optimized for minimum bandwidth. Data 
compression element 2935 supplies a formatted multi-band RTK data stream. An optional 
data link 2940, such as a radio Unk, Intemet connection, FM-subcarrier modulation, cell 
30 phone link or other transmission link, transmits raw observables of the reference receiver to 
the rover in the formatted multi-band RTK data stream. Raw observables of the reference 
receiver are processed at the rover using factorized 3+ carrier ambiguity resolution filters in 
accordance with embodiments of the invention for example as described above. 

Figure 30 illustrates an embodiment of a reference receiver processor 3005 for preparing a 
35 formatted multi-band, multi-GNSS RTK data stream, e.g., as indicated in Figure 28. Carrier- 
phase and pseudorange measurements for each of bands 1,2, . . . « of a first GNSS are 
collected, e.g., in memory blocks 3010, 3015, . , . 3020. Carrier-phase and pseudorange 
measurements for each of bands 1, 2, . . . « of a second GNSS are collected, e.g., in memory 
blocks 3025, 3030, . . . 3035. The reference receiver clock supplies a time signal to a block 
40 3040 which provides time time tags for each epoch of measurement data, e.g., in intervals of 
0.1 second. A data set per epoch is suppUed by memory blocks 3010, 3015, ... 3020 for the 
first GNSS and by memory blocks 3025, 3030, . . . 3035 for the second GNSS to an optional 
data compression element 3050. 
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Also supplied to data compression element 3050 is a known or calculated position 3045 of the 
reference receiver, which is used in RTK positioning to calculate a vector from the reference 
receiver to the rover. Data compression enables transmission overhead (time and/or cost) to 
be minimized, or enables the data rate to be optimized for minimum bandwidth. Data 
5 compression element 3050 suppUes a formatted multi-band RTK data stream for the first 
GNSS and the second GNSS. Data compression element 3050 supplies a formatted multi- 
band RTK data stream. An optional radio hnk, Intemet connection, FM-subcarrier 
modulation, cell phone link or other transmission link, transmits raw observables of the 
reference receiver to the rover in the formatted multi-band RTK data stream. Raw 
10 observables of the reference receiver are processed at the rover using factorized 3+ carrier 
ambiguity resolution filters in accordance with embodiments of the invention for example as 
described above. 

Figure 31 shows a mode of operation in accordance with embodiments of the invention in 
which a single reference station is employed. A reference receiver 3105 supplies reference- 

15 station data 3110, for example a formatted multi-band RTK data stream as provided by 
processor 2905 in Figure 29 or processor 2035 in Figure 30. Reference-station data are 
provided for example via a suitable data link. A rover receiver 3115 comprises suitable 
elements for receiving reference-station data 3110, such as a data link antenna 3120 and 
datalink signal reception and demodulation electronics 3125, and supplying reference-station 

20 data 3 1 1 0 to a rover processor 3130. If reference-station data 3 1 05 are compressed for 
efficient transmission via the data link, a data decompression element 3125 within rover 
receiver 31 15 decompresses the reference-station data for use by one or more fiuther 
processing elements within processor 3130. 

Rover receiver 31 15 further includes an antenna 3140 for receiving GNSS signals having 
25 three or more carrier frequencies and suitable signal reception and demodulation electronics 
3 145 for producing rover-reciever 3+ carrier phase and pseudorange data 3 145 for use in 
rover processor 3130. Rover processor 3130 includes a multi-band carrier-phase and 
pseudorange synchronization element 3150 which synchronizes rover-receiver data 3145 with 
reference-station data 31 10 for use in a multi-band RTK position computation element 3155. 
30 A computation element 3 1 60 uses reference-station data 3 1 1 0 to compute the reference- 
station position and provide this with reference-station description to position computation 
element 3155. An optional computation element 3165 computes atmospheric and satellite 
corrections for used by position computation element 3155. Computation element 3155 
computes RTK position fixes of rover receiver 3115 and supplies as output data RTK position 
35 and status information 3 1 70. 

Figure 32 shows a mode of operation in accordance with embodiments of the invention in 
which a plurality of reference stations is employed in a network. Reference receivers 3205, 
3210, 3215 each supply reference-station data, for example a formatted multi-band RTK data 
stream as provided by processor 2905 in Figure 29 or processor 2035 in Figure 30, to a 

40 network RTK server 3220. Network RTK server 3220 generates satellite and atmospheric 
corrections and supplies these with the reference-station data for example as a stream of 
network RTK multi-and corrections 3225 for use by one or more rover receivers such as rover 
receiver 31 15. Network RTK multi-band corrections 3225 are provided for example for use 
in a virtual-reference-station mode or in a broadcast mode. The virtual-reference-station 

45 mode is comparable to the mode of Figure 31 as the reference-station data is combined by 
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network RTK server 3220 to produce a set of data simulating a reference station (e.g., 
simulating single reference-station data 3110 from reference receiver 3105) for a declared 
location such as the location of rover receiver 3115. The broadcast mode supplies reference- 
station data from each of reference receivers 3205, 3210, 3215 for processing in rover 
5 processor 3130. In either mode, corrections 3225 optionally include error models such as 
ionosphere and troposphere models calculated from the reference-station data. 

Figure 33 shows a mode of operation in accordance with embodiments of the invention in 
which a plurality of individual reference stations is employed. Reference receivers 3305, 
3315, 3325 supply respective reference station data 3310, 3320, 3330 to a rover receiver 

10 3300. The data are provided for example via a suitable data link. Rover receiver 3300 

includes suitable elements for receiving the reference-station data, such as a data link antenna 
3335 and datalink signal reception and demodulation electronics 3340, and supplying the 
reference-station data to a rover processor 3345. If the reference-station data are compressed 
for efficient transmission via the data link, a data decompression element 3350 within rover 

15 receiver 3345 decompresses the reference- station data for use by one or more fixrther 
processing elements within processor 3345. 

Rover receiver 3300 further includes an antenna 3370 for receiving GNSS signals having 
three or more carrier frequencies and suitable signal reception and demodulation electronics 
3375 for producing rover-reciever 3+ carrier phase and pseudorange data 3380 for use in 

20 rover processor 3345. Rover processor 3345 includes a respective multi-band carrier-phase 
and pseudorange synchronization and reference information element 3355, 3360, 3365 for 
each of the reference stations. Elements 3355, 3360, 3365 synchronize rover-receiver data 
3380 with reference-station data for use in a multi-band RTK position computation element 
3390. Computation element 3390 computes RTK position fixes of rover receiver 3300 and 

25 supplies as output data RTK position and status information 3395. 

Those of ordinary skill in the art will reaUze that the following detailed description of the 
present invention is illustrative only and is not intended to be in any way limiting. Other 
embodiments of the present invention will readily suggest themselves to such skilled persons 
having the benefit of this disclosure. Reference will now be made in detail to 
30 implementations of the present invention as illustrated in the accompanying drawings. The 
same reference indicators will be used throughout the drawings and the following detailed 
description to refer to the same or Uke parts. 

In the interest of clarity, not all of the routine features of the implementations described herein 
are shown and described. It will be appreciated that in the development of any such actual 

35 implementation, nvmierous implementation-specific decisions must be made to achieve the 
developer's specific goals, such as compliance with application- and business-related 
constraints, and that these specific goals will vary from one implementation to another and 
from one developer to another. Moreover, it will be appreciated that such a development 
effort might be complex and time-consuming, but would nevertheless be a routine 

40 undertaking of engineering for those of ordinary skill in the art having the benefit of this 
disclosure. 
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In accordance with embodiments of the present invention, the components, process steps 
and/or data structures may be implemented using various types of operating systems (OS), 
computer platforms, firmware, computer programs, computer languages and/or general- 
purpose machines. The methods can be run as a programmed process running on processing 
circuitry. The processing circuitry can take the form of numerous combinations of processors 
and operating systems, or a stand-alone device. The processes can be implemented as 
instructions executed by such hardware, by hardware alone, or by any combination thereof. 
The software may be stored on a program storage device readable by a machine. 
Computational elements, such as filters and banks of filters, can be readily implemented using 
an object-oriented programming language such that each required filter is instantiated as 
needed. 

Those of skill in the art will recognize that devices of a less general-purpose nature, such as 
hardwired devices, field programmable logic devices (FPLDs), including field programmable 
gate arrays (FPGAs) and complex programmable logic devices (CPLDs), apphcation specific 
integrated circuits (ASICs), or the like, may also be used without departing firom the scope 
and spirit of the inventive concepts disclosed herein. 

In accordance with an embodiment of the present invention, the methods may be implemented 
on a data processing computer such as a personal computer, workstation computer, mainfi-ame 
computer, or high-performance server running an OS such as Microsoft® Windows® XP and 
Windows® 2000, available from Microsoft Corporation of Redmond, Washington, or 
Solaris® available fi^om Sun Microsystems, Inc. of Santa Clara, California, or various versons 
of the Unix operating system such as Linux available firom a number of vendors. The 
methods may also be implemented on a multiple-processor system, or in a computing 
environment including various peripherals such as input devices, output devices, displays, 
pointing devices, memories, storage devices, media interfaces for transferring data to and 
firom the processor(s), and the like. Such a computer system or computing environment may 
be networked locally, or over the Intemet. 

Methods and apparatus in accordance with embodiments of the invention include but are not 
limited to: 

1 . A method of processing a set of GNSS signal data derived from signals having at least three 
carriers, comprising: 

a. Applying to the set of GNSS signal data a geometry filter using a geometry carrier- 
phase combination to obtain an array of ambiguity estimates for the geometry carrier- 
phase combination and associated statistical information; 

b. Applying to the set of GNSS signal data a bank of ionosphere filters using a 
geometry-free ionosphere carrier-phase combination to obtain an array of ambiguity 
estimates for the ionosphere carrier-phase combination and associated statistical 
information; 

c. Applying to the set of GNSS signal data at least one bank of Quintessence filters 
using a geometry-free and ionosphere-free carrier-phase combination to obtain an 



Attorney Docket A-1403 



page 55 



array of ambiguity estimates for the geometry-free and ionosphere-free carrier-phase 
combination and associated statistical information; 

d. Applying to the set of GNSS signal data at least one code filter using a plxirality of 
geometry-free code-carrier combinations to obtain an array of ambiguity estimates for 
the code-carrier combinations and associated statistical information; and 

e. Combining the arrays of a., b., c. and d. to obtain a combined array of ambiguity 
estimates for all carrier phase observations and associated statistical information. 

2. The method of 1 , ftirther comprising computing a user position from the combined array. 

3. The method of 1 or 2, wherein the set of GNSS signal data comprises GNSS signal data of a 
single epoch. 

4. The method of any of 1-3, fiirther comprising computing a user position from the combined 
array with a combination of float ambiguities and fixed ambiguities by applying an integer- 
least-squares procedure and a validation procedure. 

5. The method of any of 1-3, further comprising computing a user position from the combined 
array with float ambiguities, and computing a user position from the combined array with a 
combination of float ambiguities and fixed ambiguities by applying an integer-least-squares 
procedure and a validation procedure. 

6. The method of any of 1-5, wherein applying to the set of GNSS signal data at least one code 
filter comprises applying to the set of GNSS signal data at least one code filter using a 
plurality of geometry-free and ionosphere-free code-carrier combinations to obtain an array of 
ambiguity estimates for the code-carrier combinations and associated statistical information. 

7. The method of any of 1-5, wherein applying to the set of GNSS signal data at least one code 
filter comprises applying to the set of GNSS signal data one code filter for each carrier using a 
plurality of geometry-free and ionosphere-free code-carrier combinations to obtain an array of 
ambiguity estimates for the code-carrier combinations and associated statistical information. 

8. The method of 7, wherein the code filters are mutually orthogonal. 

9. The method of any one of 1-8, wherein the set of GNSS signal data is derived from signals 
received from transmitters of a first navigation system, the method further comprising 
processing a second set of GNSS signal data derived from signals having at least two carriers 
and received from transmitters of a second navigation system by: 

f. Applying to the second set of GNSS signal data a geometry filter using a geometry 
carrier-phase combination to obtain an array of ambiguity estimates for the geometry 
carrier phase combination and associated statistical information; 

g. Applying to the second set of GNSS signal data a bank of ionosphere filters using a 
geometry-free ionosphere carrier-phase combination to obtain an array of ambiguity 
estimates for the ionosphere carrier-phase combination and associated statistical 
information; 
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h. Applying to the second set of GNSS signal data at least one bank of Quintessence 
filters using a geometry-free carrier phase combination to obtain an array of ambiguity 
estimates for said carrier phase combination and associated statistical information; 

i. Applying to the second set of GNSS signal data at least one code filter using a 

5 plurality of geometry-free and ionosphere-free code-carrier combinations to obtain an 

array of ambiguity estimates for the code-carrier combinations and associated 
statistical information; and 

wherein d. further comprises combining the arrays of f., g., h. and i. with the arrays of a., 
b., c. and d. to obtain a combined array of ambiguity estimates for all carrier phase 
10 observations and associated statistical information. 

10. The method of 9, wherein a. and f. are carried out by applying the first-mentioned set of 

GNSS signal data and the second set of GNSS signal data to a single geometry filter to obtain 
a single array of ambiguity estimates for the geometry carrier-phase combination and 
associated statistical information. 
15 11. The method of 9 or 10, wherein the first navigation system has a first number of carrier 

frequencies and observables and the second navigation system has a second number of carrier 
frequencies and observables which differs from the first number of carrier frequencies and 
observables. 

12. The method of any of 9-1 1, wherein applying to the second set of GNSS signal data at least 
20 one code filter comprises applying to the second set of GNSS signal data at least one code 

filter using a plurality of geometry-free and ionosphere-free code-carrier combinations to 
obtain an array of ambiguity estimates for the code-carrier combinations and associated 
statistical information. 

13. The method of any of 9-1 1, wherein applying to the set of GNSS signal data at least one code 
25 filter comprises applying to the set of GNSS signal data one code filter for each carrier using a 

plurality of geometry-free and ionosphere-free code-carrier combinations to obtain an array of 
ambiguity estimates for the code-carrier combinations and associated statistical information. 

14. The method of 13, wherein the code filters of h. are mutually orthogonal. 

15. The method of any of 1-14, wherein the GNSS signal data comprises data collected at a 
30 plurality of reference receivers. 

16. The method of any of 1-15, wherein the GNSS signal data comprises data collected at a 
plurality of roving receivers. 

17. The method of any of 1-16, wherein the GNSS signal data comprises data generated from a 
network of reference stations. 

35 18. The method of 1 7, wherein the data generated from a network of reference stations comprises 

at least one error model. 
19. Apparatus for carrying out the method of any of 1-1 8. 
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20. Apparatus for processing a set of GNSS signal data derived from signals having at least three 
carriers, comprising: 

a. A geometry filter using a geometry carrier-phase combination to obtain from the set 
of GNSS signal data an array of ambiguity estimates for the geometry carrier-phase 
combination and associated statistical information; 

b. A bank of ionosphere filters using a geometry-free ionosphere carrier-phase 
combination to obtain from the set of GNSS signal data an array of ambiguity 
estimates for the ionosphere carrier-phase combination and associated statistical 
information; 

c. At least one bank of Quintessence filters using a geometry-free and ionosphere-free 
carrier-phase combination to obtain from the set of GNSS signal data an array of 
ambiguity estimates for the geometry-free and ionosphere-free carrier-phase 
combination and associated statistical information; 

d. At least one code filter using a plurality of geometry-free code-carrier combinations to 
obtain from the set of GNSS data an array of ambiguity estimates for the code-carrier 
combinations and associated statistical information; and 

e. A combiner to produce from the arrays obtained by the filters of a., b., c. and d. a 
combined array of ambiguity estimates for all carrier phase observations and 
associated statistical information. 

21. The apparatus of 20, further comprising a position-computation element to compute a user 
position from the combined array. 

22. The apparatus of 20 or 21, wherein the set of GNSS signal data comprises GNSS signal data 
of a single epoch. 

23. The apparatus of any of 20-22, further comprising a position-computation element to compute 
a user position from the combined array with a combination of float ambiguities and fixed 
ambiguities by applying an integer-least-squares procedure and a validation procedure. 

24. The apparatus of any of 20-22, further comprising a position-computation element to compute 
a user position from the combined array with float ambiguities, and computing a user position 
from the combined array with a combination of float ambiguities and fixed ambiguities by 
applying an integer-least-squares procedure and a validation procedure. 

25. The apparatus of any of 20-24, wherein said at least one code filter uses a plurality of 
geometry-free and ionosphere-free code-carrier combinations to obtain an array of ambiguity 
estimates for the code-carrier combinations and associated statistical information. 

26. The apparatus of any of 20-24, wherein said at least one code filter uses a plurality of 
geometry-free and ionosphere-free code-carrier combinations to obtain an array of ambiguity 
estimates for the code-carrier combinations and associated statistical information. 

27. The apparatus of 26, wherein the code filters are mutually orthogonal. 
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28. The apparatus of any of 20-27, wherein the set of GNSS signal data is derived from signals 
received from transmitters of a first navigation system, the apparatus further comprising the 
following elements for processing a second set of GNSS signal data derived from signals 
having at least two carriers and received from transmitters of a second navigation system: 

e. A geometry filter using a geometry carrier-phase combination to obtain from the 
second set of GNSS signal data an array of ambiguity estimates for the geometry 
carrier phase combination and associated statistical information; 

f. A bank of ionosphere filters using a geometry-free ionosphere carrier-phase 
combination to obtain from the second set of GNSS signal data an array of ambiguity 
estimates for the ionosphere carrier-phase combination and associated statistical 
information; 

g. At least one bank of Quintessence filters using a geometry-free carrier phase 
combination to obtain from the second set of GNSS signal data an array of ambiguity 
estimates for said carrier phase combination and associated statistical information; 

h. At least one code filter using a plurality of geometry-free and ionosphere-free code- 
carrier combinations to obtain from the second set of GNSS signal data an array of 
ambiguity estimates for the code-carrier combinations and associated statistical 
information; and 

wherein said combiner produces from the arrays obtained by the filters of e., f., g. and h. 
and the arrays obtained by the filters of a., b., and c. a combined array of ambiguity 
estimates for all carrier phase observations and associated statistical information. 

29. The apparatus of any of 20-27, wherein the set of GNSS signal data is derived from signals 
received from transmitters of a first navigation system, the apparatus further comprising the 
following elements for processing a second set of GNSS signal data derived from signals 
having at least two carriers and received from transmitters of a second navigation system: 

f. A bank of ionosphere filters using a geometry-free ionosphere carrier-phase 
combination to obtain from the second set of GNSS signal data an array of ambiguity 
estimates for the ionosphere carrier-phase combination and associated statistical 
information; 

g. At least one bank of Quintessence filters using a geometry-free carrier phase 
combination to obtain from the second set of GNSS signal data an array of ambiguity 
estimates for said carrier phase combination and associated statistical information; 

h. At least one code filter using a plurality of geometry-free and ionosphere-free code- 
carrier combinations to obtain from the second set of GNSS signal data an array of 
ambiguity estimates for the code-carrier combinations and associated statistical 
information; and 
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wherein the geometry filter uses a geometry carrier-phase combination to obtain from the 
first-mentioned set of GNSS signal data and from the second set of GNSS signal data an 
array of ambiguity estimates for the geometry carrier phase combination and associated 
statistical information, and wherein said combiner produces from the arrays obtained by 
the filters of f., g. and h. and the arrays obtained by the filters of a., b., c. and d, a 
combined array of ambiguity estimates for all carrier phase observations and associated 
statistical information. 

30. The apparatus of 28 or 29, wherein the first navigation system has a first number of carrier 
frequencies and observables and the second navigation system has a second number of carrier 
frequencies and observables which differs from the first number of carrier frequencies and 
observables. 

3 1 . The apparatus of any of 28-30, wherein said at least one code filter of h. uses a plurality of 
geometry-free and ionosphere-free code-carrier combinations to obtain from the second set of 
GNSS signal data an array of ambiguity estimates for the code-carrier combinations and 
associated statistical information. 

32. The apparatus of any of 28-3 1, wherein said at least one code filter of h. uses a plurality of 
geometry-free and ionosphere-free code-carrier combinations to obtain from the second set of 
GNSS signal data an array of ambiguity estimates for tiie code-carrier combinations and 
associated statistical information. 

33. The apparatus of 32, wherein the code filters of h. are mutually orthogonal. 

34. The apparatus of any of 1-33, wherein the GNSS signal data comprises data collected at a 
plurality of reference receivers. 

35. The apparatus of any of 1-34, wherein the GNSS signal data comprises data collected at a 
plurality of roving receivers. 

36. The apparatus of any of 1-25, wherein the GNSS signal data comprises data generated from a 
network of reference stations. 

37. The apparatus of claim 36, wherein the data generated from a network of reference stations 
comprises at least one error model. 

While embodiments and applications of the present invention have been shown and described, it will 
be apparent to those skilled in the art having the benefit of this disclosure that many more 
modifications than mentioned above are possible without departing from the inventive concepts 
herein. The invention, therefore, is not to be restricted except in the spirit of the appended claims. 
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